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Bemerkungen zur Stabilisierung eines physikalischen Pendels. I 
Von Wolfhart Haacke in Braunschweig 


Es werden die Stabilisierungsmöglichkeiten eines im Aufhängepunkt vertikal periodisch erschülterten 
Pendels gesucht. Es wird gezeigt, daß bei Betrachtung des Pendels als starren Körper nur die normale 
Gleichgewichtslage und die um an verschobene Lage stabil sein können, während bei Berücksichtigung der 

ı  Blastizität sich eine schräge stabilisierte Lage ergibt. 


The author investigates the possibilities of stabilizing a pendulum that is periodically disturbed by shocks 
acting vertically at the point of suspension. If the pendulum is considered a rigid body, only Ihe usual position 
of equilibrium and the opposite position prove stable, while, taking account of the elasticity, an oblique stable 
position is resulting. 


Les possibilites sont discutees de stabiliser un pendule pousse dans son point de suspension verticalement 
et periodiquement. L’auteur montre, que le pendule considere comme corps rigide n’a aucune position d’equi- 
libre que sa position normale et celle deplacte par n; en temant compte de l’ölasticite on oblieni une position 
oblique. 


B padore HCCHEAYIOTCH BOBMOKHOCTU CTAÖHAH3BANHH MAATHHKA, MOABEPTAIIIETOCH B TOYKE 
HOABeCa BEPTHKAJIBHBIM TIEPHOAUYECKUM COTPACEHHAM. Ecnu mastuur P&CCcMATPHBATb KaAR 
TBEPA0e TEeJIO, TO yCTOÜYHBBbIM MO3KeT ÖBITb TOABRO HOPMAJIBHO® HOJIOSKEHHE PABHOBECHA 
HJIH IIOJIOSKEHHE, AHAMETPAJIBHO IPOTHBONOJIOSKHOE TEPBOMY. IIpu y4eTe yIpyTOCTH MAATHHKY 
yCTOHYHBBIM ORA3bIBAeTCH HAKJIOHHOE TOJIOSKEHHE. 


1. Aufgabenstellung 


Wenn ein physikalisches Pendel in seinem Aufhängepunkt vertikal periodisch erschüttert 
wird, so ergeben sich zusätzliche Gleichgewichtslagen. Diese Erscheinungen sind in mehreren 
Arbeiten untersucht worden: [1] bis [9]!). Im wesentlichen wollen wir an die Arbeiten von 
Erdelyi [2] und Klotter [4] anknüpfen. 

Nach Klotter ergibt sich für eine gewisse (kritische) Erschütterungsfrequenz ejne Stabili- 
sierung des Pendels in. der oberen — sonst labilen — Lage, die mit wachsender Erschütterungs- 


kritische Frequenz, doch folgert er auf Grund» 
seiner Theorie, daß sich mit wachsender Erschüt- 
terungsfrequenz die Gleichgewichtslage immer 
mehr der horizontalen Lage nähert. Es soll im 
folgenden unter anderen gezeigt werden, daß 
diese Aussagen Spezialfälle der allgemeinen Er- 
scheinungen beschreiben. 


Wir betrachten das in Bild 1 skizzierte 
Pendel. Seine Bewegungsgleichung lautet 


Bild 1 


®" + © cos ee) sin®—=0 (1). 


4 Dabei ist / die reduzierte Pendellänge, «= 2 rent die unabhängige Veränderliche, t die Zeit, 


‘g die Erdbeschleunigung und n = 4 cos 9, nt die Erschütterungsfunktion. 


le 7, Notwendige Bedingung für Gleiehgewichtslagen in den oberen Quadranten 
Wir wollen zunächst nach der Existenz kleiner (beschränkter, aber nicht notwendig peri- 


© 5 odischer) ‚Schwingungen um eine Mittellage & fragen. Diese Mittellage werde durch 


ti 


t>o 


[m ee 2.4.0) 
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definiert. Notwendige Bedingung für die Stabilisierung des Pendels ist, daß die linearisierte 
Bewegungsgleichung für kleine Schwingungen beschränkte Lösungen aufweist. Mit P®=&+9 


bei eI<Z ergibt sich aus (1): 
Be 2 Me BGE 
p +4 eos u— ern a anal 7 08 W) Sina ee (3). 


Es läßt sich sofort die Partikularlösung der inhomogenen Differentialgleichung (3) gp = — tg & 
angeben, d.h. also, es ist vr eine Konstante, und zwar für den zunächst zu betrachtenden oberen 


Halbkreis 0 <lal <ö 5 —. Nach dieser Bemerkung kann die Untersuchung auf die homogene 


Gleichung 
j gcosa 4A Rt 
(0) (808%, 4% ac0s: cos u)g=0 . NEAR N 
x . e i nl: g9c05 & A 
beschränkt werden, eine Mathieusche Gleichung, die wir durch — ana ze T c0S%=Yy 
in die Normalform @"+ (A+y cos u)p= 0 überführen. Die Beschränktheit der Lösungen 


dieser Gleichung in Abhängigkeit von A,y ist ausführlich untersucht worden [10], S. 40. Wir 
können sie nach Strutt graphisch (Bild 2) in einer A,y-Ebene darstellen. Für unsere Frage- 
stellung ist cos «>0, also A<O und y>0. Wir haben daher die Verhältnisse im zweiten Qua- 
dranten der A,y-Ebene zu untersuchen. Die einfach-schraffierten Bereiche enthalten Parameter, 
denen zwei unabhängige beschränkte Lösungen zugeordnet sind. Die Grenzkurven „enthalten“ 
eine periodische und eine unbeschränkte Lösung, 
sie gehören also noch zum instabilen Bereich. Die 
analytische Darstellung der Grenzkurven entnehmen 
wir aus Jahnke-Emde [11]. Nach Umrechnung 
der Parameter erhalten wir für 


- 7 29 
09: Ay = — 5 He 99 Pr Tan 
1 RL y° y° y 
Ge: A Are p) Be 8 4 39 384 F 
2 3 4 (4). 
1 I Yy 


R Fe ER 
Ber A EST an ar 


5 
8 Ey: Age, — 17 Parc 
Zunächst wollen wir uns mit Hilfe des Bildes 2 einen groben Überblick verschaffen. Wenn wir 
bei festem A die Erschütterungsfrequenz des Pendels erhöhen, so bewegen wir uns in der A,y- 
Ebene auf einer Geraden y = konstant von links nach rechts. Zunächst befinden wir uns im 
instabilen Bereich, bis wir auf die Kurve ce, stoßen. Wir gelangen nur dann nochmals in ein 
instabiles Gebiet, wenn diese Gerade oberhalb P die y-Achse schneidet. Den Punkt P wollen 
s 2 

wir bestimmen. Es ist dort ungefähr z + nn —0, also : cos a—=y = 0,449. Selbst für 
«—=0 muß A—=0,449 1 sein, d.h. die gesamte Erschütterungsstrecke 2A ist fast so lang wie 
die reduzierte Pendellänge, während für «>0 A noch größer sein müßte, damit die Gerade 
im Punkte P die y-Achse schneidet. 


Uns interessiert die Abhängigkeit der Lösungen von 1/n? und cos «. Daher bilden wir die 


4 
Parameterebene durch <= 7-45, 2 = und y= T cos & auf die z,y-Ebene ab. Diese 
Abbildung ist für cos&=+0 eindeutig. Dann ergibt sich aus (4): 
7. 29 1 Ya, 
= 3.1, 2° RP BEI I I I 
Deeo —=2y— str fr Kse, — Yy 2A 8 2006 
_ IT re ara a, 
Dee, y r2+5 start Ya Et Sg 
Zur graphischen Darstellung in Bild 3 benutzen wir aber an Stelle dieser Näherungsreihen die 
von Ince berechneten Tafeln [12]. Mit den dortigen Bezeichnungen ist = — 2 und y=44. 
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Zur Stabilität unseres Pendels ist 2>0 und y>0 notwendig, weiter müssen in Frage kommende 
Punkte unterhalb der Horizontalen y = 7 liegen. Wir beschränken uns jetzt auf“ < 0,44, so 
daß unsere Horizontale die y-Achse unterhalb von ce, schneidet. Die Verhältnisse für große A/l 
kann man ebenfalls leicht überblicken. Sie sollen im weiteren nicht mehr betrachtet werden. 


Bild 3 Bild 4 
) - Der schraffierte Bereich in Bild 4 stellt stabile Verhältnisse für ein festes A dar. Der Punkt 
1 
PR Fe bestimmt die kritische Frequenz n,, bei der das Stabilisierungsphänomen für a—=0 
Pat Zu einem Punkt 2,< x, d.h. n,>n,, erhält man ein Stabilisierungsintervall. Es 
ergeben sich also kleine Schwingungen um alle Lagen x mit 0 <a <a, — are cos en N. be- 
stimmt sich in erster Näherung aus 9 = 2" T = 45 ‚ also n, = BT Dies ist auch der von 


Klotter und Erdelyi gefundene Wert. Erdelyi art nach einer periodischen Lösung und 
erhält daher ce,, also die untere Grenzkurve unseres Bereichs. Diese Kurve gehört als Grenz- 
noch zum Instabilitätsbereich; stabile Lagen ergeben sich aber beliebig nahe oberhalb ce,. 
Klotter [4] fragt von vornherein nur nach Stabilität für «—=0 und erhielt daher die Gerade 


y= = für 9>22>0. Aus Bild 3 ersehen wir, daß diese Gerade zum Stabilitätsbereich ge- 


hört. Es ist aber nach unseren bisherigen Ergebnissen auch jede Lage zwischen der Klotter- 
schen Geraden und der Erde&lyischen Kurve möglich, so daß wir von einem indifferenten Bereich 
sprechen müßten. Allerdings haben wir dabei bisher die Reibung- und die Elastizität des Pendels 
sowie die Bedeutung der Mittellage und den Einfluß der Anfangsbedingungen noch nicht in 
unsere Betrachtungen einbezogen. Das soll später nachgeholt werden, vorher aber wollen wir 
noch einiges sagen über die 


3. Gleichgewichtslagen in den unteren Quadranten 
Es si Ben — a, Y=n—d=ß-+y. Statt (1) ergibt sich für Y 


N BE 7 A 11 
een 0 11a dh Eee) 
und entsprechend 
Ben —% cos Beos u) y=0 en er ee. (DR, 
Tun: 1 
geosß ae: 
& "Also ist a or 37:1 Boy Ce Pu 
Wir müssen jetzt die Verhältnisse im vierten Quadranten der /,y-Ebene untersuchen. Da 
k die „Struttsche Karte‘ zur A-Achse symmetrisch ist, können wir uns auf den ersten Quadranten 


beschränken und dabei mit (4) und Bild 2 arbeiten. Uns interessieren jetzt die doppelt-schraf- 
41 { 
fierten Bereiche des Bildes 2. Da A>0 und y>0 ist, müssen wir wegen & se und y=y in 


_ der @,y-Ebene den zweiten Quadranten betrachten (Bild 3). Man sieht leicht ein, daß für alle 
n 11* 
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A A 
», se,(v>0) die negative z-Achse Asymptote ist. 7. schneidet also je sl 


für hinreichend kleines n alle diese en, so daß — wie bekannt — Stabilitätsbereiche ; : 
mit*Instabilitätsbereichen abwechseln. Auch ° 
hier (Bild 5) wollen wir nur den Fall 


Jede Gerade = 


ce 


7 <0,44 diskutieren. Dann wäre nach die- 


ser notwendigen Bedingung für @J< || = 
das ganze Intervall 0<Sß< Z stabil. Ist 
aber 2,<%;, so en sich Stabilität für 
0<P,< B<Z . Die natürliche Gleich- 


gewichtslage ist nicht mehr stabil. Analog. 
wirken sich dann die Kurven höherer Ord- 
Be nung aus. Doch für diese Erscheinungen 
ae ist n schon recht klein, so daß die‘ 
experimentelle Bedeutung geringer ist. Für 

sehr Born n ist jedenfalls die untere Lage, wie auch Klotter [4] et stabil. 


Ein 4. Die Auswirkung der Reibung 
N Bei Berücksichtigung einer geschwindigkeitsproportionalen Reibung gilt statt (1) 
k i = 
ii - " ’ 4 g 
DB" +2xB +7 00su— Im no. : 
Analog ergibt sich für (3): 
n ‚(4 9 9 4 DEN, 
107 Harp HG 0004) 00509 (oT eosw)sina. : 
Weiter ist pp =—tga. Die homogene Gleichung lautet also: ERS 


" 9% 
pt 2x + (dma- Anen * )eosa:p=0. 


Mit 9=e-'“%:n(u) erhalten wir 


" cos & A 5 
$ ah (r ee) + Zeosa cosuln=0. 
Im oberen Quadranten mit cos «>0 ist A—= (# er je nn Er = 4 2 E 


Wir befinden uns also wieder im zweiten N. der A,y- -Ebene. 2 ‚benutzen wir 


‚die Abbildung 2 = — LRBN ‚ y=y. Dann erhalten wir wieder x — cos Ba 


4A 
man a : 


Ps nr 


also die alte Abhängigkeit. Aus 


RE N 
6 3033? 
ergibt sich jetzt 
j 4x2 7 
N TER 2y—— 
Dre 0 Yy Kr 8 
und analog 
ER) 
m +24 
4-1 
= 24 
Yy 
er a 
Va TR y2 Te 33 
2) Dabei ist 2% der Reibungskoeffizient. | 
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Das Bild 3 ändert sich also etwas: Alle Kurven werden um einen kleinen Betrag gehoben. ce, 
schneidet die %-Achse etwa in 


= VE yı-yı-7 02 jar. 


er 


Doch diese Aussagen gelten für 7. Wegen = e- “%n(w) wird beiderseits der Bereich erwei- 
tert, im Prinzip aber hiernach die Erscheinung des ‚‚indifferenten Stabilitätsbereiches‘“ durch 
die Reibung nicht beeinflußt. 


5. Die möglichen Mittellagen beim starren Körper 
Wir haben in (2) als Mittellagen 


t>o 


t 
ee lim — I B(rT)dr 
0 


definiert. Diese Definition ist mit 
. t N 
Zt 
Hu | o(T)dr—=0 


2 

L gleichwertig. Als Lösung von (3) ergibt sich nach dem Floquetschen Theorem [10], S. 57: 
} = —tga + et fu)+netef(— u). 
£ 

; 

| 

5 

} 

E 


Dabei sind c, und c, die Integrationskonstanten, u ist der charakteristische Exponent, der für 
beschränkte Lösungen rein imaginär sein muß. Daher setzen wir mit reellem » u=i». Es ist 


Ju) = f(u-+2r). Dann ist 
P=—tga-+ c, eiru 5 a, eltu oe -iru > I 


=—o r.=770 


=—tga+ 5 eiru [cq) le GALGENT, 


=—o 
Also muß 
es 
2 wem | > eiru [e, a, era ,a_„e-iv’u] du 
7 u Kiegheg! x=—-0 
B- . - , . . N « ® H 
x sein. Da wir Ya; <C LEN können, ist Summation und Integration vertauschbar. Wir 


erkennen, daß nur für ganzzahlige » eine von Null verschiedene linke Seite sich ergibt. Denn 
dann wäre 


hr 


u>o U 


„at 65) I] Sraul- te 


Hiernach wäre also für ganzzahlige » die Lage innerhalk des indifferenten Bereichs nur von den 
Be Anfangsbedingungen abhängig. Doch dieser Schluß ist nicht richtig. Ist u=ni (n ganz), so 
ist  (u+2r) = e"irp(u)=g(u). Eine Lösung ist also periodisch, die zweite daher nicht 
beschränkt. Der Fall «—=n: muß also für unsere Betrachtungen ausgeschlossen werden. Dann 
ist, falls die Lösung von (3) beschränkt ist, immer tg x == 0. Wir erhalten nur die von Klotter 
angegebene Lage. 
6. Die Elastizität des Pendelkörpers 
Wir wollen nun den Einfluß der Elastizität im Spezialfall des mathematischen Pendels 
_ untersuchen®). Dazu betrachten wir einen durch eine Federkraft biegesteif an den Aufhänge- 
_ * punkt gebundenen Massenpunkt m. R sei die Länge des ruhenden Pen- 
 — _ dels ohne Einfluß der Schwerkraft mg, d.h. bei entspannter Feder 
- (Bild 6). Ist r die,augenblickliche Länge des Pendels und ® der Winkel 
gemäß Bild 1, so beträgt zunächst beim Fehlen einer Erschütterung die 
5 potentielle Energie dieses mathematischen Pendes U— 5 (R—r)? 
+ mgrcos ®. Dabei ist o die Federkonstante. Die Werte für r und ®, 
_ für die das Potential U(r,®) zum Minimum wird, sind die Koordinaten 
“  Tiefpunkte von ‚‚Potentialmulden“, in denen Stabilität vorliegt®). 
Ku rt 


Bild 6 


U,=eo(r—R)-H-mgcos d=0, Us=—mgrsin ö=0, 


FAN 3) Di g zu dieser Untersuchung verdanke ich Herrn Prof. Dr. R. Weyrich. 
5! Oi .) Der Stab de ist in der Elastizitätstheorie noch nicht vollständig untersucht. Wir wollen 
er einfashe, ‚Definition Ye deren Hulescknaen mit den Experimenten im Einklang stehen. 


PErLz re 


“ , 
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und die Diskriminante 
A= TUT, Ur -Up=— mgorcos D— m’ sin? ®d. 
Wir erkennen drei Lösungen: 


ER AZUF r= BL Ei A, = —mge| zegeg: 

e ß 
Dun: = R+T, A,— mge(R+"2) >0. 
) Bares”. nl); A,= oR— m’g?20. 


Falls ein Extremum vorliegt, so ist dieses wegen U„> 0 sicher ein Minimum. Im Fall a kann 
für hinreichend kleines o (wenig Ähnlichkeit mit einem starren Körper) die Lage D—=0 ein 
Minimum der potentiellen Energie zur Folge 
haben, also stabil sein. Im Fallbliegt immer ein 
Minimum’ vor, die untere (natürliche) Lage ist 
auch bei Berücksichtigung der Elastizität stabil. 
Der Fall c, das Analogon zur Erdelyischen : 
Lage, ergibt hingegen für hinreichend großes o x 
ein Minimum (Ausartungsfall). 

Nach dieser Einführung soll das erschüt- 
terte Pendel betrachtet werden. Während sich 
beim starren Pendel die Bewegungsgleichungen- 
auf ein bewegtes Koordinatensystem mit seinem 
Ursprung in dem erschütterten Aufhängepunkt 
bezogen, wollen wir der Untersuchung beim 
elastischen mathematischen Pendel ein mit der Erde starr verbundenes Ruhesystem (zunächst 
kartesische Koordinaten x, y) zugrunde legen (Bild 7). Das bewegliche Attraktionszentrum 
sei 0 (0, Acos2 nt). Dann ist die potentielle Energie von m 


U= = fa: + (y— Acos2rnt”? —2R Yx? + (y— 4A cos 2rnt)?] + mgy-- const. 


i 
E 
F 


Bild 7 


VOOWEEROTANET < zu 


Zunächst soll gezeigt werden, daß wir einen kleinen Fehler begehen, wenn wir von U nur 
den zeitlichen Mittelwert betrachten. Dazu entwickeln wir U in eine Fourierreihe 


a oo 
= 1 N.a,c0os2unvt, 
er 


A Bu 


in der die a, Funktionen des Ortes sind. Wir wollen zeigen, daß die Lösung bis auf o(h) 


nur von, dem’Mittelwert von U abhängt. Dabei soll die Masse m kleine gern in der 


0a 

x, y-Ebene um ein eigentliches Minimum (z,, %) durchführen. Dann ist —— un =(.- Wir 
führen für 2—= x, und y= y, folgende Abkürzungen ein: 67 9 

1 %a 1.9 1 9a 1 da 1 da 

b _— we = m — 9 ee Se ei d er al 

GR WER Je 1a = By 2m oray” ° 2m oy? UT NO. er a °y 


bei v»>0. Dann Bit b11 des — bi, >0, b5,>0 und 5,,>0. Wegen der kleinen Schwingungen 
können wir |©— x, und |y—yo| als hinreichend klein annehmen, so daß bei NOrDachTPSSiEgzE 
höherer Glieder die Bewegungsgleichungen für m io1gpnds Form haben: 


64 bula— a) + bualy— Yo) + 2 0,005 2anni=0, ® 


Y-+ bz,(®— Ro) + bag (u U) + 2 av cos 2nnvi—0. 
Zur weiteren Ra dieses EEE von Differentialgleichungen setzen wir zur Abkürzung: 
BER bis 
A (4) = ar Das ge) = =Am— (db +5) AR (br dr — bi,), 
und hiermit ist 


oo 1 ? 
i=m+E+ > Aoan) [d12 ©2, —(b2g An? n? v?) c,,] cos Arınvt, 
v=1 N - \ 


t 


jr © 5 

j / 1 . 
FG =YWtn+ — — [da1 4,» — (bj, — 4 n? v2) c5,] cos 2unvt. 
N i Yy Yo 2 Bam)l 21 "19 ( 11 IT ) av]. e%, 
J L 
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&,n sei die Lösung der homogenen Gleichungen. Hierin sind die Reihen gleich o(5,) ‚ können 
N 
also für hinreichend großes n vernachlässigt werden. Für & und n ergibt sich dann einfach: 


Et budtben-d, ntrbattrban—t, 

ein Differentialgleichungssystem mit konstanten Koeffizienten, dessen charakteristische Glei- 
chung A(A)=0 lautet. Die Lösung hängt also tatsächlich nur von a, ab. A(A)— 0 hat nur 
reelle Wurzeln. Daher ergibt sich für m bis auf O(1/n?) eine quasielastische schwingende Bewegung 
um den Punkt x,, %,, die, da die Wurzeln von A(A)=0 im allgemeinen inkommensurabel sind, 
fastperiodischen Charakter hat. Übrigens wird durch die Tatsache, daß die in x und y auftreten- 
den zusätzlichen Fourierreihen ihrer Kleinheit wegen vernachlässigt werden können, nach- 
träglich gerechtfertigt, daß in der Entwicklung nach (x«—x,) und (y—Y,) diese in a, bis zu den 
Gliedern zweiter Ordnung, in den a, (v > 1) aber nur bis zu.den Gliedern erster Ordnung vor- 
getrieben wurde. Damit ist unsere Behauptung bewiesen, daß wir den Fehler O(1/n?) begehen, 
wenn wir von U nur den zeitlichen Mittelwert @|2 betrachten. 

Zur rechnerischen Durchführung wollen wir Polarkoordinaten verwenden und erhalten 
nach Bild 7 


= 5 (r— RP + mggcos ®; 


® wird dabei ‚„‚raumfest‘‘ gemessen, nicht aber r. Daher ist in U im raumfesten ®,g9-System 
r?— A? cos?2rnt + E —2Agcos2nntcos ® zu setzen. Führen wir jetzt die zeitliche Mit- 
) telung durch, so erhalten wir 


2x 
ee 9 Rt mggcos +8 Re 4 SE VPFEOST Ian D cos rdr 
N) 


SERIE gr ae = 
9, (R Drmgges® +, 4 1 in 8/0 2)’ 
Für A=0 ergibt sich der oben betrachtete Fall. Es ist 
2 
Us = —sin S(mag-+ Ze cos ) =); 


4° 3 4? En 
U,=mgcos ®B+o(g—R)+ ER de D—(,; N) mer; oRsin?®, 
2 n2 3 2 R 
A = —mggo cos D— m? g?sin? D — Zr cos2 D-+ en sin? Dcos ® 
4 2 P2 
— |. sin? ® 
Die drei möglichen Fälle sind jetzt folgende: 
f mg A? Ro? 
a Da; = R——, A, = m? 9 — mqaRo — = ——. 
sign im A, = — 1. Für hinreichend große o (‚fast starrer Körper‘) liegt ein Sattelpunkt vor. 


Die obere Lage ist also im Sinne unserer Definition instabil. 
a : Mg A?Ro? 
b Da, Er a: A, = m? 9’ + mg Ro — — . ——. 
Auch hier liegt für hinreichend große o ein Sattelpunkt vor, doch bleibt das Minimum länger 
erhalten als im ersten Fall. 


c) D=arc cos e Be 


Wir sehen, daß mit wachsendem o ®gegen x]2 strebt. (Hier also nicht mit wachsendem n)! 
Aus U,=0 folgt: 5 
\ { 2 m: g? q? A?o 
R 5 old R)— AroR 42 
h: Für sehr große o gilt: 

5 1 4 R 

N @?— R?+—-=0. 

Ä Wir setzen q=R-e, dann gilt 

@+2Re + Rre+ ER 


= 0% 


RT NE 


.. 


. 


» 


a h IRRESN 


& 
" 
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oder angenähert 


— R+ YR?—2 4 AR 3R+YR—-224 42 


a ex BL a8 
ERREN. A en 


Are ARE q? 112 mA gi gt 
= Tim 2m?g 1) 3 MR 


Wegen ge Rund A< g ist für hinreichend große o sicher sign A, = + 1. Weiter ist 
| a 34% ee) ( 1 za 
Im-elTamrsmt mir 


Diese ‚„‚schräge“ Lage ist also stabil im Sinne unserer Definition’). 
Wir erkennen aus dieser Betrachtung, daß die elastischen Eigenschaften des Körpers 
Anlaß sind, daß das Pendel eine andere als die Klottersche Lage einnimmt®). 


7. Zusammenfassung 


Die verschiedenen Untersuchungen über die Stabilisierung eines physikalischen Pendels 
mit periodisch vertikal erschüttertem Aufhängepunkt konnten bisher noch nicht in Einklang 


miteinander gebracht werden. Die hier im wesentlichen benutzte Methode wurde zuerst von 
Klotter [4] angewandt, doch nur bei der speziellen Fragestellung, wann das Pendel für «=0 
stabilisiert wird. C 


Durch die Berücksichtigung der verschiedenen Einflüsse konnten die unterschiedlichen 
bisher beobachteten Erscheinungen erklärt werden. Durch verschiedenartige Versuchsbedin- 
gungen und besonders durch unterschiedliche Bauweise der Pendel wurde bisher sowohl die 
schräge als auch die vertikale Mittellage experimentell bestätigt. In einer kleinen Versuchsreihe 
beobachtete der Verfasser beide Lagen. Eine schräge Lage ergab sich, sobald sich kleine Biege- 
schwingungen des Pendelkörpers einstellten. Durch sorgfältig ausgewählte Anfangsbedingungen 
konnten diese Biegeschwingungen ausgeschaltet werden. Dann wurden die Aussagen der Klotter- 
schen Theorie bestätigt. Dabei wurde die Lage der Erd&lyischen Theorie in dem Sinne labil, 
daß von ihr aus das Pendel sowohl nach unten als auch nach oben ‚‚fallen‘‘ konnte. Die tech- 
nische Schwierigkeit dieser Experimente liegt darin, daß bei den notwendigen hohen Erschütte- 
rungsfrequenzen ein Spiel der Gelenke und elastische ROHWEIE UNGEN nur schwer ausgeschaltet 
werden können. 


Die Lösung der Differentialgleichung (3) ergab sowohl für den oberen als auch für den 
unteren Quadranten bei hypothetischer Einführung der Mittellagen (2) indifferente Parameter- 
bereiche für stabile Lösungen. Doch konnte mit Hilfe des Grenzwertes (2) gezeigt werden, daß 
nur die Lagen «=0 und «= in Frage kommen. Bei dieser Untersuchung, die durch Reibungs- 
effekte im Prinzip nicht geändert wird, ist die linearisierte Differentialgleichung als Aus- 
gangspunkt gewählt worden. Weiter ist hierbei die Elastizität des Körpers nicht berücksich- 
tigt. Bei den Experimenten beobachtet man außerordentlich kleine Schwingungen um die 
jeweilige Mittellage, sog. Zitterschwingungen. Daher brauchen wir auf den Unterschied der Lö- 
sungen (1) und (3) nicht einzugehen. 4 


‚Neu ist die Betrachtung des elastischen Einflusses. Dieser läßt gerade eine schräge Lage 
“als ausgezeichnet erkennen. Je nach der Größe dieser elastischen Einflüsse wird die schräge 
oder die Klottersche vertikale Lage eingenommen werden. 


Die Untersuchung dieser Arbeit beschränkt sich aber auf vertikale lineare Oszillation des 
Aufhängepunktes, während mehrere der Arbeiten [1] bis [9] allgemeinere ebene Oszillationen 
zulassen. Schon bei der Betrachtung einer linearen nicht vertikalen Oszillation ergibt sich die 
Erschwerung, daß sich 9, nicht mehr auf eine Konstante reduzieren läßt. Doch mit Hilfe des 
Satzes, daß es bei nichtperiodischen Lösungen von (3)n immer eine periodische Partikularlösung 
von (3) gibt, sind analoge Schlüsse auch im Abschnitt 5 durchführbar. Es ergeben sich aber im 
allgemeinen andere Gleichgewichtslagen. — Herrn Prof. Dr. Weyrich möchte ich für An- 
regungen und Gespräche danken, durch die er das Entstehen dieser Arbeit gefördert hatıa = gs 


x % a 


Pr 


5) Es muß aber darauf hingewiesen werden, daß bei dieser Theorie die Frequenzabhängigkeit TER zum 
Ausdruck kommt. Daher kann aus dieser Betrachtung nicht geschlossen werden, daß sich das Pendel mit 
wachsender Erschütterungsfrequenz der horizontalen Lage nähert. Auch in den später zu er 
Experimenten konnte dies nicht beobachtet werden. 

6) Die Spezialisierung auf das mathematische Pendel ist schon deshalb eine Einschrär g, weil ier 
nur Längsschwingungen, aber keine Biegeschwingungen auftreten können. Die analoge Behan Re des | 
kalischen Pendels dürfte sehr schwierig sein, so daß wir uns auf Se Untersuchung E ag 

d 


das Wesentliche klar erkennen läßt. 


re ca " 4 


ie 
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Graphische Integrationsverfahren für die Meißnerschen 


) - Linienbilder bei Differentialgleichungen 1. Ordnung 
Fi. Von Rudolf Ludwig in Braunschweig 


Während Meißner mit seinem -Verfahren eine Differentialgleichung 1. Ordnung nicht direkt inte- 
grieren konnte, werden hier neue Methoden angegeben, mit denen dies möglich ist. Außerdem kann auf 
einfache Weise eine Fehlerabschätzung durchgeführt werden, die sonst bei graphischen Verfahren meist 
unterlassen wird. ; 


While the method of Meissner is not applicable to differential equations of the first order, the author 
gives new methods that are filling this gap. Furthermore, the estimate of the error is easily accomplished that 
is generally ommitted in ihe case of graphical methods. 


La methode de Meissner ne permettant pas l’integration directe d’une Equation differentielle du 
premier ordre l’auteur donne des procedes nouveauz pour resoudre ce problöme. De plus une estimation 
d’erreur, d’ordinaire ommise dans le cas d’un procede graphique, ici se fait simplement. 

Tax kax meTon MmeücHepa HeUPMHMeHUM MIA UHTErPuUPOoBaHHuA AubepeHnmasbHbIX YpaBHeHuH 
HePBOTO MOANKA, TO ABTOP IpemIaracT HOBbIe METONBI, SANOTHAMINHE ITOT ıpodes. Kpome 
TOTO YKASBBAETCH IIPOCTOÜ CIOCOO OMECHKH IIOTPEINHOCTH, YTO OÖBIYHO NPH TpahmyeckoM 
HHTerPHPOBAaHHN HE TelaeTcH. 


Einleitung 

In einer früheren Arbeit!) wurden für bestimmte Integrale Methoden angegeben, um mit 
den Meißnerschen Linienbildern graphisch gute Näherungswerte zu bestimmen. Bei den 
 Meißnerschen Verfahren ist die Integration von DGln.?) zweiter und höherer Ordnung durch- 
] führbar, eine DGI. 1. Ordnung konnte nicht direkt behandelt werden. Nachdem Doetsch [1a] 
und Schulz[1b] für DGln. 1. Ordnung ganz andere Wege gegangen waren, sind hier— aufbauend 
auf den Methoden für bestimmte Integrale — Konstruktionen mit einem rasch konvergierenden 
‘ Iterationsverfahren entwickelt worden. Eine Fehlerabschätzung geht von einer von Tollmien 
[2a], [2b] stammenden Methode aus, die in eine für verschiedene Anwendungen besonders ein- 
fache Form gebracht wurde. In einer späteren Arbeit sollen die hier dargelegten Verfahren auf 

DGin. beliebiger Ordnung ausgedehnt werden. 


1. Allgemeines zur näherungsweisen Integration und Restabschätzung von 
Differentialgleicehungen 1. Ordnung 


P- m f 2 

: Es soll die explizite DGl. 1. Ordnung 

Es: | 20.2 ll) 
_ behandelt werden. Gesucht wird die Lösung, die der Anfangsbedingung 

Be a nalen) 


1) „Analytische Untersuchungen und konstruktive Erweiterungen zu den graphischen Integrations- 
hren von Meißnerund Grammel.“ Z. angew. Math. Mech. 31 (1951) 120—130, Hinweise auf das 
i n und die Formeln dieser Arbeit werden durch eine vorgesetzte „1“ gekennzeichnet. 

DGlI(n) = Differentialgleichung(en), Gl(n) = Gleichung(en). 


De 
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genügt. Von der Funktion f(x, y) werde vorausgesetzt, daß sie 1. beschränkt sei: 


f(x, y)|&.&: 22: ea EN (2) 


in einem Gebiet &: |e — z,|<a, |y—y,|<b und 2. der folgenden Lipschitz-Bedingung 
genüge: 


Kaya Fa ym)IsMyı pl... ner. er 8). 


Diese Bedingung ist sicher dann erfüllt, wenn man die weitere Voraussetzung macht, daß die 
partielle Ableitung der Funktion f nach y existiert und daß diese beschränkt sei 


of * 
B ER EN (3%) 
Die formale Lösung von (1) lautet: 
yazutlflar@)de. #2. 0 BE ar 2 (4). 


Praktisch kann man die Gl. (4) durch ein Iterationsverfahren lösen. Bezeichnet man die Anzahl 
der Iterationsschritte durch oberen Index, so ergibt sich 


yo) — y, +/1% yo (a) dr were Ar eeE a Eee (>: 


Dieses Picardsche Iterationsverfahren konvergiert, wenn (2) und (3) erfüllt sind. 


Zunächst werde an eine Abschätzungsformel erinnert (s. z.B. [3]), mit deren Hilfe 
man den Unterschied zwischen Näherungslösung und exakter " Lösung erfassen kann. Ersetzt 


man in (1) die Funktion f(z, y) durch eine Näherungsfunktion f(z, y), so entsteht eine neue DGl., 
die der Anfangsbedingung 


It) = Io yo = ö 
genügen soll. Sie habe die Lösung y(x), also ist 


I a HA-FEIH) H+HR IE). 20: (6) 


Diese Lösung (x) soll zur Ausgangsfunktion einer Iterationsfolge der gegebenen DGI. (1) ge- 
nommen werden. Es stehen sich also die beiden GlIn. gegenüber: 


9) 10= [Fa a) = F1@ Hm)dr + [9a s(2))dz SE Wh} 


* 
Ey =). Mm v8) dw ee re ee BE ee (8). 
To - 
Setzt man das Iterationsverfahren bei der letzten Gl. fort, so findet man 
Bu 
N) — ya yrTDla))da.-.- m... (9). 
%e v - » 


Aus (7) und (8) folgt: 
> yD (a) — J(8)= — 6 — [ dla, VEN) BR EEE 10). 
%o 
Die Iterationsfolge (9) konvergiert nach den Voraussetzungen (2) und (3)nach y(x). Das Integral 
der Gl. (10) kann man abschätzen: - 
Be % ’ 
h. Pe s@)dz <o a N ER ae 1) 
N 
Also: Fu > 
UELIORZIOJEGEIE C Sr ee 2. (109 
Aus der Differenz der beiden Gln., die aus (9) für n und n—1 entstehen, findet man die Ab- 
schätzung: 
\ T 
|y® (2) — yrd @)| <M [|yr® (x) — y"%(a)| de 
) % - 


\ 


Z.angew, Math. Mech. 


A „Z. angew. Math. Mech 


N ee Th a NE 
re I Ba f h 
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Durch rekursives Einsetzen ergibt sich die auch von Tollmien [2] angegebene Abschätzungs- 
formel: 


y—y|<jöjerie ir em). A WAL. (18% 


Gilt für die Näherungsformel Y(&) dieselbe Anfangsbedingung wie für die DGI. (1), ist also 
lv, So:ist: 


- y—1< Gern 1) NEE dan): 


Für die weitere Fehlerabschätzung eines Näherungsverfahrens soll der Methode von Tollmien 
gefolgt werden, aber sie soll in eine etwas einfacher zu handhabende Form gebracht werden. 


Die exakte Lösung (4) der DGl. wird durch eine Näherungslösung y(x) ersetzt, indem das 
auf der rechten Seite stehende Integral näherungsweise bestimmt wird. Man kann dazu solche 
Näherungsverfahren, wie sie in I entwickelt Pe etwa in folgender Form angeben: 


Int mh > anf (Anıaı Ins) BEN a 1A): 


Eine solche Gl. läßt sich im allgemeinen an ein numerisches oder graphisches Iterations- 
verfahren lösen. Ist jedoch &,—0, so kann man die Gl. direkt zur Bestimmung von Yy„--z aus 
den vorhandenen Werten benutzen. Die Lösung der Näherungs-DGl. (6) kann man schreiben: 


Kr E7 var %,+kh 


7, 
2 Sf (in+ıs In+4) ER NS (15). 
Von der Näherungslösung y(x) sind die diskreten Werte y,für &,— 2, vh bekannt. Ebenso 
ist die Funktion 
fa, 9) = F(a,) 


nur an diesen Stellen bekannt. Aus obiger Gl. folgt: 
%&,rkh %ytkh 
% 32 H@))da—h af ann); 7 I(®, 9(2)dx) 
%,trkh 
3 = lan) — 1 F(a) de. 


Der letzte Ausdruck ist aber das Restglied der Integration, also der Integrationsfehler 
&,tkh 


} ET RE Be N: (16) 


und daraus: 
A| 3 
O, = kh R, (2, y(x)) ’ 
wobei das Restglied mit der gefundenen Näherungslösung y(x) zu bilden ist. Nach Gl. (11) gilt: 


| R, IC} E 


Be 2 


E20 
\[9(& 90) da] < |o— | 9 < |2— zo] |Onlmaz = 
%o 
oder 
1 
ir 9= 7, |Rulmaz RER E) N RR BF (17). 
_ Für die DGI. 1. Ordnung gilt somit folgende Fehlerabschätzung einer Näherungslösung y: 


ie Imas (rien — 1) ng RN GEN 


1 


y»—7|s 


Hieraus ist sofort zu ersehen, daß sich der Fehler aus zwei Faktoren zusammensetzt, nämlich 
‚dem Integrationsfehler [4] für einen Schritt und der Fehlerfortpflanzung durch den gesamten 
Integrationsbereich. 
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2. Das Tangentenverfahren ohne Iteration (Bild 1) 
Um die DGl. (1) mit der Anfangsbedingung (1*) graphisch zu integrieren, ind also zu- 
nächst nur je ein Punkt der Polaren O und der ersten Polaren CO, und damit auch nur ein Punkt 
des Linienbildes Z bekannt. Wenn die Polare der ersten Ableitung bekannt wäre, brauchte man 
nur eins der Integrationsverfahren I (2,2 bis 2,4) an- 
zuwenden. Nimmt man zunächst einmal an, es seien 
nach Festlegung der Schrittweite % die Punkte P, und 
P; für i=0,1,..., k schon bekannt — die Kon- 
struktion solcher Ausgangspunkte soll in Abschn. 5 
besprochen werden —, so läßt sich sofort das Tangenten- 
verfahren 1(2,3) zur Integration über ein Intervall2 Ah 
anwenden. Nach Gl. I(4) ist folglich 


Yr+ı = Ir-ı + 2, f, sin e 
oder, wenn man die Schrittweite so wählt, daß Be 
e=arcsinh, | 8 
Bild1. Tangentenverfahren für DGln. 1. Ordnung ehr Je 92 hf (& 9). : pe 


Die Lösung erfolgt also in diesem Falle ohne Iteration. Mit dem soeben konstruierten Wert 57.11 
berechnet man 


Ir = f (+1 Ir+1) 
und wiederholt die Konstruktion nun mit den Punkten P,, P;+ı und P; und Prıı. 


Zur Durchführung der Restabschätzung findet man mit der zweiten Gl. I (13) Zür das 
Intervall 24: 


EN RE re 
o- m | 6 If (Ey) EL (19). 
und damit 


Verse NR HM - 2 a 


Dabei ist der erforderliche Anfangswert y, als fehlerlos vorausgesetzt, ist dieser jedoch mit dem R, 
Fehler ö behaftet, so gilt: = 


ET ET BER Eee 


Iy—J| Ss del el u —_ es a ea EBEN BE Re 
M 
h 3. Das Sehnenverfahren mit Iteration re 
$ Br. Für das Sehnenverfahren angewandt auf die DGI. (1) ergibt sich nach Gl. I (7) aus Aue 1 u 
; (2,2): : a 
RR Ir = I + Lf(& 9) + f&rrı» I) 5 T ee 
7, wobei — arc 85 h gesetzt wurde. Hat man außer den Punkten P, und 14 noch einen ersten. 


2 
Näherungswert von aERN: so kann man daraus 


er Ha 3) Ge art, 


berechnen und damit 93 P}')ı zeichnen. Dies bedeutet, daß man obige Gl. durch Trade : ” 
lösen .kann, indem man stets die Sehnen im Linienbild zieht. Die Bedingung dafür, daßdas 
Verfahren konvergiert, ist bekanntlich, wenn die zu lösende Gl. die Form IF, Fly) Bat 


h Ei 


Br; 


>. RD 
< 2. angew. Math. Mech. 


X er ae ver ” u Er FDA | 
A EN FAR Be ee yr ey! er 
vr Tr ee « au y. 


REN E ? 
al > Er a i f RT 


RR | 
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Dadurch ist also eine Bedingung für die Schrittweite festgesetzt. Um im ganzen Integrations- 
bereich gleichgünstige Verhältnisse der Konvergenz zu haben, muß man beachten, daß das 
Produkt AM bei stark veränderlichem M ungefähr konstant gehalten wird. Hierdurch wird 
nach Collatz eine ‚natürliche Schrittweite‘ [5]festgelegt. Beieiner ‚Schrittkennzahl‘“ MA in 
der Größenordnung von etwa 0,2 wird man erfahrungsgemäß meist eine so günstige Iterations- 
folge bekommen, daß schon der zweite Schritt graphisch kaum in Erscheinung tritt. Um einen 
ersten Näherungswert zu erhalten, verwendet man das im letzten Abschnitt angegebene 
Tangentenverfahren. 


Für die Restabschätzung ergibt sich aus der ersten Gl. I (13) der Integrationsfehler: 
1 (T h? em 2 
9= ,|Rn = Da (x, (El) Be le ra N (22), 
und daraus die Abschätzungsformel: 


Be ES A (2, I(&)) Imaz a 7 a PERS BEN (23). 


Der Fehler kann also gegenüber dem Tangentenverfahren durch Anwendung des Sehnen- 
verfahrens mit Iteration nur auf die Hälfte vermindert werden. 


4. Das kombinierte Sehnen-Tangenten-Verfahren mit Iteration (Bild 2) 


Wie in den vorhergehenden zwei Abschnitten seien wieder die drei Punkte P,, P; und Q; 
für}=0,1,..., k bereits bekannt, wobei die Schrittweite A festgelegt worden sei} Nach dem 
Tangentenverfahren (Abschn. 2) bestimmt man einen nullten Näherungswert pa, und be- 
rechnet mit der DGl. (1) einen Näherungswert von pP: Die drei Punkte Pi_1, Pr; Pat 
sind Punkte der ersten Polaren der gesuchten Integralkurve. Integriert man diese erste Polare 
mit P;_ı als Anfangspunkt über das Inter- 
vall2% nach dem Sehnen-Tangenten-Ver- 
fahren I(2,4), so erhält man einen ersten 
verbesserten Wert pi) T, Wieder berechnet 
man mit der DGl. (1) den verbesserten Wert 
P®, und wiederholt das Verfahren solange, 
bis innerhalb der Zeichengenauigkeit der 
Punkt Pf), sich nicht mehr ändert. In 
Bild 2 'ist bereits PP, = P%,ı. Nach. 
1(2,4) gilt für das Iterationsverfahren 
Gl. 1(14) 


>= Je er (f(x, Ir) 


+f(&+1,9xH1))tg h+ 45 (9, 9) sin h]. 
Aus dieser Gl. folgt als Bedingung für die 
Konvergenz 


= | J tg h < I, Bild 2. Sehnen-Tangenten-Verfahren für DGIn. 1. Ordnung 
oder Mtgh<3. 
.Da h genügend klein ist, kann man auch einfacher dafür schreiben: 
Ri | | en Eee NEU 
Aus der Gl. I(15) für den Integrationsfehler folgt sofort: 
ER RNIT Im OD) 


h h- - Die Fehlerabschätzung für das Sehnen-Tangenten-Verfahren lautet also 
er . 


IHREN HF E NH IE De (6). 
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Dieser Fehler ist im Vergleich zu dem des Sehnenverfahrens (23) größenordnungsmäßig um den 
Faktor Ah? kleiner, außerdem besagt die Konvergenzbedingung, daß eine größere Schrittweite 
zulässig ist. Bei vorgegebener Fehlerschranke muß man also bei dem Sehnenverfahren das 
Intervall in eine größere Anzahl Schritte unterteilen als bei dem nur wenig mehr Zeichenarbeit 
erfordernden Sehnen-Tangenten-Verfahren. 

Bei der Anwendung des allgemeinen Sehnen-Tangenten-Verfahrens über ein %-fach unter- 
teiltes Intervall kh nach 1 (2,5) müssen also % Punktepaare P,, ..., Paz _ı und Pr, - - , Par 
bekannt sein. Nach der Tangentenkonstruktion Gl. I (17) gewinnt man wie vorher einen 
nullten Näherungswert P}),; und berechnet mit der DGl.(1) 9%, um den Punkt P,), 
zeichnen zu können. Nun läßt sich auch der Sehnenzug nach Gl. I (16) Banane! und durch 
Bildung des Mittelwertes gemäß Gl. I (18) erhält man den verbesserten Wert PN) . und wieder- 
holt nun die Konstruktion solange, bis sich der Wert von Y„.„innerhalb der Zeichengenauigkeit 
nicht mehr ändert. Aus den Gln. I(16), (17) und (18) folgt als Konvergenzbedingung für das 
Iterationsverfahren: 


5k—6 |df . h <1 
(BER) Fa 
oder 
h a 
Mtg SEHE N IE EEE (27), 
wofür man bei kleinem h auch schreiben kann: 
12 (k—1) 
MI< DR RR BER (27*) 
Man erhält für k—=2 wie in Gl. (24) Mh<3 und im -Grenzfall für k— oo 
MRS ED NE ER DER ER (28). 


Diese letzte Bedingung kann man, da die rechte Seite am kleinsten ist, auch unabhängig von k 
stets zur Festsetzung der Schrittweite fordern. 


5. Die Bestimmung von Ausgangswerten (Bild 3) 


Durch das Anfangswertproblem der DGI. (1) und (1*) sind von den Polaren C und ©, und 
dem Linienbild Z der gesuchten Lösung nur je ein Punkt bekannt. Zur Durchführung der in den 
vorhergehenden Abschnitten dargelegten Verfahren benötigt man aber wenigstens je zwei Punkte. 
Der Satz2aus 1(2,1) legt die Tangentenrichtung der Polaren in sehr einfacher Weise fest. Für 

i . benachbarte Punkte kann man also einen Näherungswert 
auf der Tangente annehmen. Dadurch ist also auch ein 
Näherungswert von y’ berechenbar. Mit diesen zwei Punk- 
ten läßt sich das Sehnenverfahren mit Iteration anwenden. 
Man muß nun zwei Fälle unterscheiden: 

a) Will man die Lösung ausschließlich mit dem Seh- 
nenverfahren durchführen, so genügt also die Festlegung 

Bild 3. Bestimmung von Ausgangswerten des Ausgangspunktes P, auf der Tangente der Polaren. 

Nur, wenn für diesen Punkt bei der gewählten Schritt- 

weite die Konvergenz der Iteration zu schlecht wird, muß man zur Unterteilung (besonders 
Halbierung) der Schrittweite übergehen. 

b) Will man das Sehnen-Tangenten-Verfahren anwenden, so muß man beachten, daß die 
Genauigkeit dieses Verfahrens größer ist als die des Sehnenverfahrens. Daher muß man durch 
Halbierung oder evtl. noch weitere Unterteilung den Fehler des nach der Sehneniteration wie 
unter a) bestimmten Punktes möglichst herabdrücken, denn der Fehler des ungenauen Anfangs- 
punktes pflanzt sich durch das ganze Verfahren fort, wie der erste Summand der Gl. (13) zeigt. 

Der Konstruktion der Ausgangswerte muß man die erforderliche Sorgfalt angedeihen 
lassen, da durch die Fehlerfortpflanzung der Fehleranteil de ==“ Jeicht den eigentlichen Ver- 
fahrensfehler übertreffen kann. | 

Außer der hier dargelegten Methode, einen ersten Näherungswert auf der Tangente der 
Polaren festzulegen, kann man bei dem Iterationsverfahren auch von einem geschätzten 
Werte ausgehen. Dies ist gerade am Anfang allerdings etwas unsicher; leichter ist dieses 
Schätzen, wenn man schon ein Stück der Kurve kennt, also nur eine Schrittweite extrapoliert. 
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6. Zusammenfassung für das Sehnen -Tangenten -Verfahren 


\ DGI.: v=f(%y) 
Anlangswerte: yo= ya)» = ao Yo) 
Integrationsbereich: N 
Festlegung der Schrittweite Ah, unter Beachtung 
der Konvergenzbedingung: Mh<3. 


1. Schritt: Bei halber Schrittweite: Pin auf der Tangente der Polaren in P, kon- 
struiert. 


Berechnet: 9% aus der DGI. 


2. Schritt: PgP/% nach dem Sehnenverfahren integriert: P\',. 
Berechnet: DAR) aus der DGlI. 
Wiederholung des 2. Schrittes bis das Iterationsverfahren innerhalb der 
Zeichengenauigkeit konvergiert. 


3. Schritt: Bei ganzer Schrittweite: nach dem Tangentenverfahren P\” konstruiert. 


Berechnet: z0 aus der DGl. 


4. Schritt: Nach dem Sehnenverfahren P, Pi» pP,” integriert. 
Durch den Mittelwert p® bestimmt. 
Berechnet: y) aus der DGI. s 
Wiederholung des 4. Schrittes (Iterationsverfahren). 

5. Schritt: Wie 3. Schritt für das folgende Intervall x, bis z,: ergibt P)). 


6. Schritt: Wie 4. Schritt: PP, PP, ..., P, 
— usw. — 
7. Beispiele 
7,1. Allgemeines 


Für eine vorgelegteDGl. (1) mit den Anfangswerten (1*) und dem geforderten Integrations- 
bereich «< 2 <b muß man sich unter Beachtung der Konvergenzbedingung für das Iterations- 


verfahren für eine Schrittweite h = er entscheiden, wobei n die Anzahl der Schritte ist. 


Nach der Größe der Zeichnung muß man dann noch den Maßstab wählen, in dem y und y’ 
aufgetragen wird. Da man x als Winkel im Bogenmaß aufzutragen hat, empfiehlt es sich, wenn 
man mehrere solche Zeichnungen anzufertigen hat, auf Transparentpapier ein Strahlenbüschel 
mit einem Winkel zwisehen den Strahlen von 0,05 im Bogenmaß aufzuzeichnen. Zu große oder 
zu kleine Winkel kann man durch affine Transformation der Variablen x der DGl. vermeiden, 
also z.B. £= cz (s. Beispiele). 

Eine prinzipielle Schwierigkeit bei der Fehlerabschätzung besteht darin, daß man höhere 
Ableitungen z. B. f!Y(x,y) der gegebenen Funktion f(x, y) braucht. Obwohl das Differenzieren 
an und für sich keine Schwierigkeiten macht, ist doch der Rechenaufwand groß und die Formeln 
werden umständlich. Ein anderer gangbarer Weg ist der, daß man die Differenzenquotienten 
als Näherungswerte benützt, dabei muß man allerdings in Kauf nehmen, daß die höheren Diffe- 
renzen stark streuen. Meist wird man allerdings graphisch den Verlauf der Differenzen glätten 
können (evtl. Wiederholung des Verfahrens bei der nächsten Differenz). 


7,2. Beispiel für das Tangentenverfahren ohne Iteration 


Es soll die DGI. | Ye y? 
mit der Anfangsbedingung: (el 
für den Bereich << 0,25 


integriert werden. Vorteilhaft ist es, die Variable « affin zu transformieren (Vergrößerung des 
Winkelbereiches): 
Durch &=nx erhält man 


day .1 [8° n 
2=.(atr) Is‘<s;- 


R ne BRPR oe re u De > 


RER \ Te ee N 
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Es wurde bei der Originalzeichnung gewählt (Verkleinerung auf etwa Y/,): RR TR ER % 
1. die Schrittweite i = 5 Hana 0,0125; . 
2. die Zeicheneinheit e = 200 mm. = 


Die Konstruktion lieferte folgende Werte, wobei die Berechnung der Funktionswerte mit einer 
Quadrattafel [6] und dem Rechenschieber durchgeführt wurde. 


n | ® y Maine Se ® | y | fa) 

0 0 1 0,3183 *) 5 0,125 ARTS TR 0,417 3 

== 0,0125 1,013 0,3268 6 0,150 1,171 0,444 

l 0,025 1,025 0,3350 FR 0,175 1,203 0,471 : 

2 0,050 1,053 0,3542 8 0,200 , 1,245 0,507 >. 

3 0,075 1,083 0,3746 9 0,225 1,283 0,540 - E 

4 0,100 1,108 0,3946 10 |; 0,250 1,327 i 0,581 & 
Der erste Wert (*) ist bei halber Schrittweite auf der Tangente der Polaren in P, konstruiert Bi 
worden. 


Beispiel 1: 
J:x?*y? yl0)-1 


Tangentenverfahren 
(ohne /teratıon) 


.--... 
mans, 


ae 4 
r rIA 
- 


Ans 


Einschließlich des Ablesefehlers kann also das Ergebnis in folgender F 
EL a a 
| Yıo == 1,33 + z gie 4 
Der Verfahrensfehler ist also in diesem Falle größer als di 
nauigkeit. 0 ı " 


> x 2 re a Be RS 
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Ben. 73. Beispiel für das Sehnen-Tangenten-Verfahren mit Iteration 
 - Für das gleiche DGl.-Problem wie im vorhergehenden Abschn. 7,2 soll das Sehnen-Tan- 


 genten-Verfahren mit Iteration in der Weise durchgeführt werden, daß die Schrittweite verdoppelt 
= 


wird, also h = 0° Man erhält folgende Wertetabelle: 
2 n | ® | y I ® | 3 |" sen) 
RE N A ETESTE EEE REEEIETT Ta 1 BE ES Sa TEE RER APRES E23 A Ser AA 
0 0 1 0,3183 2 | 0,100 1,109 0,3948 
Ri. Ze — 0,0125 1,013 0,3268 3 0,150 1,174 0,447 
= = 0,025 1,026 0,3268_ 4 0,200 1,245 0,506 
1 0,050 1,052 0,3358 5 0,250 1,341 0,593 
| Innerhalb der Zeichengenauigkeit genügt stets eine einmalige Iteration. Die Li pschitzkon- 
Ep stante lautet:  ; = - £ 
| | | 5 of 2, 
- SION man SRRTE n/A 


Berspiel2: 
Vexy3 Jo) 


Sehnen- Tangenten-Verfahren. 
(mit Iterationen) 


‚diesem Wert ist also die Konvergenzbedingung 
RN E, % iR le, FR j F AN a 
z 20 < 


a AR 4 Lat ? 7 

von 0,1 mm sind also für den letzten gefundenen Wert 9; die , 
1. Die Fehlerschranken darf man nicht zu eng wählen, da 
henungenauigk eiten hinzukommen, die sehr stark von der 


abhängig s 


ER 


12 


Er SE 5 Fa ER a a a 
8 Y t ee a - 
, Br n a er. ze ch U 
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Für die Fehlerabschätzung kann man einer Differenzentafel der 7 folgende Maximalwerte 
entnehmen: 


22 250,0257 [AR fan 0,028, flnae = 0:5: 
Daraus folgt für die Differentialquotienten: 
AR ee 
WE 7SE Bel al! If iss = 12. 


Die Formel für die Restabschätzung liefert folgende Schranken: 


7 
19% — ls apa 190.7 40 + 12 + 54] (e* — 1) = 0,0003. 


Der nichtberücksichtigte Fehler der Ausgangswerte spielt hier keine Rolle. 


8. Über ein numerisches Verfahren 
Die in I Abschn. 2,6 abgeleiteten Integrationsformeln (I, Tabelle 2) lassen sich auch sehr 
einfach zur numerischen Integration von DGln. verwenden. Sind bereits die Werte Y, und 5, 
für v»—=0, 1,..., n—1 bekannt, so findet man mit der betreffenden Formel für J; : 


yi0) — Ina + J%. 
Mit diesem Wert berechnet man nach (1) y/(%. Jetzt kann man am einfachsten J* berechnen 
und dann den gemittelten Wert J* benutzen: - | 
yU— In. + J* y 
und dann mit diesem Wert weiter iterieren. 7 } 
Ausgangswerte kann man in ähnlicher Weise wie bei den graphischen Verfahren aus der 
ersten Zeile der Tabelle I, 2 gewinnen, indem man mit J* beginnend für das halbe Intervall 
einen Näherungswert bestimmt. 
Die in diesen Formeln enthaltene Integration unter Benutzung der Simpsonschen Regel 
soll an einem von Willers [7] behandelten Beispiel gezeigt werden: 
Für die DGI. 


2 rec ee ee ee 


G 


= —Yay—1 
mit der Anfangsbedingung ’ 


‚t)=2 
erhält man: r 
k | x (0) I DR SL ee BE El, ya) yi2) 
0 1 2 = ’ n— 4 1,4 1,57381 1,57418 — 
—_ 1,05 1,95 1,94942 —_ 5 1,5 1,46420 1,46454 _ 
1 LT 1,89768 1,89774 _ usw, i 
2 1,2 1,79143 1,79188 1,79187 17 2,7 0,44416 0,44425 0,44424 
3 1,3 1,68324 1,68364 E= 18 2,8 0,40381 0,40383 — 
Striche in den Spalten 4und 5 bedeuten, daß keine Änderungen gegen die vorhergehenden Spalten eintreten. 
Die Werte sind die gleichen, wie sie von Willers angegeben wurden, aber es sind höchstens 
zwei Iterationen nötig. Die Fehlerabschätzung liefert: € 


Be IYıs — Yıs| < 1.10. 


Schrifttum . 


[la] G.Doetsch: Über die graphische Integration von Differentialgleichungen erster Ordnung. Z. angew. 
Math. Mech. 1 (1921), S. 464—466. 
[1b] W. Schulz: Über das Meißnersche Integrationsverfahren für Differentialgleichungen erster Ord- 
nung. Dtsch. Math. 6 (1941), S. 271—276. \ 
[23] W. Tollmien: Über die Fehlerabschätzung beim Ada msschen Verfahren zur Integration gewöhn- 
licher Differentialgleichungen. Z. angew. Math. Mech. 18 (1938), S. 83—90. 
[2b] A- Fricke: Über die Fehlerabschätzung des Adam sschen Verfahrens zur Integration gewöhnlicher | 
Differentialgleichungen 1. Ordnung. Z. angew. Math. Mech. 29 (1949), S. 165—178. 4 
[3] R. Rothe: Höhere Mathemavik. Teil III. 1. Auflage. $83. S. 177. Leipzig 1935. 
J. Horn: Gewöhnliche Differentialgleichungen. 3. Auflage. $11. Gl.(11). Leipzig und Berlin 1937. 
[4] W. Tollmien: a.a.0. Die Gl. (13) ist auf S. 85 als Gl, (11) richtig angegeben, aber bei der Durch- 


rechnung in $5, 8. 90 muß es in den Schlußformeln (56) und (56a) statt 77 se“ ei] heißen deM «==, 


[5] L. Collatz: Natürliche Schrittweite bei numerischer Integration von Differentialgleichungssystemen. 
Z. angew. Math. Mech. 22 (1942), S. 216—225. ker 
[6] Z.B. in der bekannten Ga ußschen 5stelligen Logarithmentafel. Wr hr ; 
[7] Fr. A. Willers: Methoden der praktischen Analysis. $33,2. S.316. (Tabelle 2, S.319.) Bein 
u, Leipzig 1928. ) | . a PS Ee 
Eingegangen: 19.6. 1950. \ 


ı 
f 
a 0 I ı u 5 


a 
y 

L 
J 


1 u Ze HEHE ze 


Z.angew, Math. Mech. 
Bd. 31 Nr. 6 Juni 1951 


Kleine Mitteilungen 179 


KLEINE MITTEILUNGEN 


Zum Gebrauch des harmonischen Analysators 
von Mader-Ott'!). 

Ich will hier zeigen, wie man 

1. durch eine bequeme, weil kongruente Umzeich- 


‚ nung der mit der Periode 2/ periodischen Kurve 


y= f(x) den vom Fahrstift des Analysators zu be- 
schreibenden Weg möglichst abkürzen, insbesondere 
das Zurückführen auf der z-Achse ersparen kann, 

2. wie man damit die Anzahl der Umfahrungen und 
der zugehörigen Planimeterablesungen, die notwendig 
sind, um nach einem von E. Nyström?) angegebenen 
Verfahren Fourierkoeffizienten zu bestimmen, für die 
keine Zahnräder vorhanden sind, und auf Perioden- 
strecken > 36cm zu analysieren, auf die Hälfte her- 
abdrücken kann, 

3. wie man dabei mit zwei zusätzlichen Umfah- 
rungen auch die dem Nyströmschen Verfahren bis- 
her nicht zugänglichen Koeffizienten a, und b, in der 
Entwicklung 


= 2 
f)-F+ 0,008 + b, sin” + ag cos —- 
23 bysin at r. 


bekommen kann. 


1. Abkürzung des Fahrstiftweges 


Aus der Theorie des Mader-Öttschen Analysators?) 
folgt unmittelbar: Sind f(x) und g(x) zwei graphisch 
gegebene Funktionen mit der gemeinsamen Periode 
21, dann erhält man mit dem Analysator die Fourier- 
koeffizienten von f(x) — g(x) genau so wie die von 
f(x), wenn man nur den Fahrstift statt auf der x- 
Achse auf der Kurve y= g(x) zurückführt. Die so 
entstehende Umfahrung wollen wir die ‚„Umfahrung 
f(x) — g(z) auf der Strecke (0,2 1)‘ nennen. Als An- 
fangs- und Endpunkt nehmen wir in der Regel einen 
Punkt der Strecke vom Punkt (0, g(0)) zum Punkt 
(0, f(0)). Die genannte Theorie zeigt auch, daß der 
Analysator für je zwei Zahlen p, 9, fürdeO <p<g< 
21 gilt, bei der Umfahrung auf der Strecke (p, g) die 
Integrale 


q x 
Ju — ga)] cos "7? da 
und ee Al), 
Al ü NIE 
Bi IT 
nl Lie) — Kajlsin TE da 


soweit die zugehörigen Zahnräder vorhanden sind. 


Bildi1 


 *ı) Ein Teil der folgenden Ausführugenn ist sinngemäß auf 
„andere Analysatoren übertragbar. 
2)E. Nyström, Über den Gebrauch des harmonischen 
- Analysators Mader-Ott. Soc. Scient. Fenn. Phys. Math. IX, 14 
(1937), referiert in den unter ?) angeführten Werken. 
») Vgl. etwa Meyer zur Capellen W., Mathematische 
Instrumente 1941, 8. 216ff. und Willers Fr. A., Mathematische 
Instrumente 1943, 8.178ff. 


Ist die Funktion f(x) in der Periodenstrecke (0,2 I) 
zu analysieren, dann zeichnen wir (Bild 1) die Funk- 
tionen f(®), f(21 — x), — f(21— x) für0<x<lein. 
Zunächst sei (bei der normalen Einrichtung des Ana- 
lysators) 2,5cm <21< 36cm. Wir stellen den Ana- 
lysator so über das Zeichenblatt, daß die Strecke 
(0,22) der Schwenkbereich des Fahrarmes ist. 


Wir schreiben die Fourierkoeffizienten 


21 
An = + [re cos 


0) 


N ITXC 
dx 
n eh 


2 
— 2 [Ifte) + f121- ajleos 


ı 
0 


NICK 
l 


dx 


21 
b, = fs) sin dx 
0 


21 
1 Dh: 
— je — f21— 2)]sin — da 


Der Analysator sei mit den Zahnrädern für n = 
1,2,3, --%k ausgestattet. Dann bekommen wir a, 
nach (1) und (1’) bei Verwendung des n-ten Zahnrades 
und „Einsatz C‘ des Planimeterstiftes mit Hilfe der 
Umfahrung f(x) + f(21— x) auf der Strecke (0, 7). 
In Bild 1 ist das die Umfahrung ADEFGHA, wir 
wollen sie kurz die Umfahrung C nennen. Ebenso be- 
kommen wir db, mit demselben Zahnrad und dem Pla- 
nimetereinsatz $S durch die Umfahrung f(x) — f(21— x) 
auf der Strecke (0,1), im Bild ADECBA, welche 
kurz die Umfahrung S heißen soll. a,ergibt sich, wenn 
man unmittelbar mit dem Planimeterstift (Planimeter- 
einsatz P) die Umfahrung C' ausführt. 

Alle diese Umfahrungen liefern übrigens dasselbe, 
wenn die Strecke (0, I) nicht gerade die linke Hälfte 
des Schwenkbereichs des Fahrarms ist, sondern irgend- 
wo in diesem Schwenkbereich weiter rechts liegt. 
Wichtig ist nur, daß der Planimeterstift beim Ein- 
satz C für x = 0',,oben‘, beim Einsatz $ ‚links‘ liegt. 

Weil Addition einer Konstanten zu f(x)nur a,/2 um 
sie vermehrt, die anderen Fourierkoeffizienten aber 
überhaupt nicht ändert, kann man durch Addition 
einer geeigneten Konstanten zu f(z), d.h. durch gün- 
stige Wahl der x-Achse den in der y-Richtung vom 
Fahrstift zurückzulegenden Weg möglichst abkürzen. 
Ist f(x) in x = 0 stetig, dann ist A = B und es wird, 
wenn man die z-Achse durch Alegt, A= B=H. 


2. Zum Nyströmschen Verfahren 


a) Für gerades n = 2 m sind 


I 
2 ) * 
20, [ve + f(21 — x)] cos — da 
0) 5 ä 
. (2) 
2 bn = lu — f(21 — 2)] sin ie 
2 


Fourierkoeffizienten der in (0, !) als f(») + f(21 — %) 
und f(z) — f(21 — x) definierten und mit der Periode I 
periodisch fortgesetzten Funktionen. Der Analysator 
liefert sie durch Analyse auf der Periodenstrecke (0, !) 
mit dem Zahnrad m. Dabei darf 2,5cm </<s 36cm 
sein, so daß man hierbei für f(x) bis zu Perioden von 
21= 72cm kommt. Hierfür bringen wir die Strecke 
(0, I) in die Mitte des Schwenkbereiches des Fahrarms. 
Diese Lage sei im folgenden immer beibehalten. 
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b) Die a„ und d„ für ungerades n = 2m — 1 be- 
kommen wir, indem wir auch in diesen Fällen die 
Integrale (1’) mit dem Nyströmschen Verfahren be- 
handeln. Wie stellen dieFahrarmlänge auf 2mI/» und 
setzen das m-te Zahnrad ein. Dann umfahren wir 
zur Bestimmung von a, nach Einsetzen des Planimeter- 
stiftes in © auf der Strecke (0, !) f(x) + 21 — 2), zur 
Bestimmung von b„ mit dem Planimeterstift in S auf 
derselben Strecke f(x) — f(21 — x). Dabei achten wir 
darauf, daß beim Beginne des Umlaufs, d. h. bei = 0, 
je nachdem wir a, oder 5b, ermitteln wollen, der Ein- 
satzpunkt des Planimeterstiftes „oben“, bzw. links 
steht. Diese Einsätze mögen Cm und Sm heißen. 
Multiplizieren wir dann noch die entstehenden 
Planimeterdifferenzen aus den von Nyström an- 
gegebenen Gründen mit m/n, dann erhalten wir a„ und 
b„n. Das geht, solange das m-te Zahnrad vorhanden 


ist und ei auf dem Fahrarm einstellbar ist. Da von 
den Zahlen „_"— für m= 2,3,.- 2/3 die größte 


ist und auf dem Fahrarm 3714 gerade noch eingestellt 
werden kann, kann man auf diese Weise noch auf 
Perioden von 56 cm analysieren, 


3. Bestimmung von‘a, und b.. 
Wir verstehen unter y(z) und ö(x) die folgenden 


H Funktionen: 
ie (2) sin 22 für OSe sa mit der Periode 
ie; Y@)T\—sin2z „na<s<e<s2nf 2n(Bild2), 
Br ö(2) = y(2 + 7/2) (Bild 3). 
(a y 
re 
Er 
E 2m > 
Br: | 
Bi IX) 
et ö 
IE Bild 2. 
N Als ihre Fourierentwicklungen berechnet man leicht 
, e 8 cos (2 m+ 1)« 
Y y(z) srl rger 
$ 2 (2m — 1)(2m +3) 
Bi e u 2-7 32T __c08 5% 
B- wÄlrS. ab ang 
F $ _ecos7a ) 
“ 5-9 ? 
“ & Ei). 
8 (— 1)” sin (2m + De 
ö(2) = — ET RTT 
: 52 (2m — 1) (2m+ 3) 
I -2[- sinz sin3z  sindx 
k a 3 RR 
sin Tx 
5:9 
Daraus gewinnen wir 
37 3cos3z  3cos5r 
car yet 175 ir a7 r 
3 cos 7x2 
5:9 
(4). 
Ta. ie Sin. Sch 3sindr 
a 1-5 SE A 
3sin 7x 
u + ... 
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Damit wird 


2 
db, + | na)sin "® a2 
F 22 . 
In TE 5) 
T sr) han 


Die Integrale | 


2ı 
1 [annfies 
f) 
ı 


2ırı 


2 
7 un + sl aylsin-TX de, 
«= [s5(%)a: 


0 

. i ; 
uw (21 216(%) er 
1/2 


- 2. [T-fe) + K21— 2)]sin 
0 


INT ; .£ 


Te 


R | 
+ ve SABHR 2)] in 22 de 


werden vom Analysator mit dem Planimetereinsatz 
S, geliefert, und zwar erhalten wir g durch die Um- 


fahrung C, s durch eine Umfahrung X, die wiram 


deutlichsten durch die Umfahrung ABCDECDA in. 
der einfachen Zeichnung Bild 1 beschreiben. 
Auf diese Weise bekommen wir a 


- 


_3n 3 eg 
ST Tee ge 


re 


5 ıl I. » 4 
—= 0,589 g9 +06, +70 +75 + Rn: 
A 
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berechnen. Die nachfolgenden Glieder zu vernach- 
lässigen ist durch die folgende Abschätzung gerecht- 
fertigt. 

Es sei 


Er Terenerra 
BALLEE ILL. 


Dann et 
2, 
Au 1 Meolgernerrn 
) cos EFMREL...|ar 
3af TE 
= [ut ar (7) da, 
wo 


a T% 
fr )=Z+ a, c08 7% 7: 


TT 


re b,sin 7 


re En 


+ Agy4g C08 


. (2k+2)nz 
+ + dDopıa ein ER FIRE 


21 
1 Prem 
kJ I 


ist, erhalten wi mit Hilfe der Schwarzschen Un- 


gleichung®) = 
21 
Dat 
} (1 As nee ae 
ö 
| 97? 
Faa9, Migte 
Ebenso ergibt sich für 
3n 3 3 IE fe 
NE Sees juni Br a Sup zu ae Fr 
z \ 3 dog+ı 
E Be REEL AN, 
E Br GT DRrrN 
| Per Be 22 
Ban.» 1 a 
(di = Br) SS 37) MD) Rrızl)] da 
8 {) 
Fr? ; gr2 
A 35 Minte: 


ee 
Die Fehler der Näherungswerte A,, und B,, von a, 
er und b, sind also absolut genommen nicht größer als 


KL 2 e a a 
ins ki, Fouriersche Reihen, $. 133. Samml. Göschen 1930, 
2 ve Pet 


I Ten N Ve 
"ule) o(e)dz. SS u)? ds [va] de; vgl. W. Rogo- 
v% re “5 
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1,67 Mggıg, WO Mgyı, der mittlere Fehler (im 
Sinne der Fehlerquadrate) von Sax +2 (7) ist. Da dieser 


nicht größer, in der Regel aber viel kleiner ist als die 
Fehler von fg;ı5 ist und man diese schon in Kauf 


nimmt, wenn man fg,,, als Ersatz für f nimmt, ist 
es gerechtfertigt Ag,,, und B,,,, an Stelle von a, 
und b, zu nehmen; z.B, ergibt s.ch bei der Funktion 


za für 0sSe<s<ıa 
k(2)= 
(22 -:) für nr ses2n, 


welche die Fourierentwicklung 


2 cos3% cosd® 
EEE nl ee La TwErT Torte 
hat, durch rechnendes Verfolgen der Definition von 
A,ı Aıı = — 0,9997, während a,= — 1 ist. Das 


Maximum des Fehlers der Näherungsfunktion k,g(2) 
wird an der Stelle x angenommen. Dort ist k(rz) — 
2,464, kj(n) = 2,423, der Fehler also ungefähr 
100 mal so groß wie der von A... 1,67 M , = 0,012. 


Übersicht 
Die Mitte des Schwenkbereiches 


1. 21<S 36 cm. 
liegt bei = 1. 


Fahrarm- 
einstellung 


Einsatz des 
Planimeter- 
fahrstiftes 


Umfahrung 


Die Planimeter- 
messung ist zu 
multiplizieren 

mit 
2. 36cm <21<S56cm. Die Mitte des Schwenk- 
bereiches liegt bei # = 1/2. 


Fahrarm- i 
einstellung 
Einsatz des 
Planimeter- S, 
fahrstiftes 
Umfahrung X 
Die Planimeter- 
messungistzu | I A. BR 
multiplizieren 1 
\ mit 
Innsbruck. L. Vietoris. 


Entfernungsmittelwerte bei der Ellipse. 

Anläßlich einer Analyse des astronomischen Grund- 
begriffes der ‚„„Mittleren Entfernung der Erde von der 
Sonne“ untersucht Emersleben!) den Mittelwert 
aller Entfernungen eines Ellipsenpunktes von einem 
Brennpunkt, :wobei sowohl über den Polarwinkel als 
auch bei konstanter Flächengeschwindigkeit über die 
Zeit gemittelt wird. 


1) O0, Emersleben: Die mittlere Entfernung Sonne-Planet, 
Astr. Nachrichten 275, Heit 5/6 (1947) S. 263, 
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Verallgemeinert kann man das kinematische Problem 
der Bestimmung der mittleren Entfernung eines mit 
vorgegebener Geschwindigkeit auf einer Ellipse um- 
laufenden Punktes vom Mittelpunkt oder einem 
Brennpunkt aus betrachten. Wir unterscheiden im 
folgenden die drei Fälle: 1. konstante Winkelgeschwin- 
digkeit, 2. konstante Flächengeschwindigkeit, 3. kon- 
stante Bahngeschwindigkeit. Es seien r und @ bzw. o 
und y Radiusvektor und Polarwinkel bezogen auf 
einen Brennpunkt bzw. auf den Mittelpunkt. Es ist 
dann 


» b 
r= ——— la e= ———— (lb), 
1+ kcosp (1m) 1 — k’cos?y 
wenn 
b? b° 
= k= 1—- — Be} 
he: ; (2) 


Halbparameter und numerische Exzentrizität be- 
zeichnen. Die Entfernungsmittelwerte sind dann bzw. 
gegeben durch 


Fahre... (3a) 
0 
4 
ERNST. 
= | etwioyaı (3b) 
oder auch‘ i 
en 
Fr Map. . (48) 
0 
2n . 
rel 
= na . (4b), 
0 


wenn & bzw. y die Winkelgeschwindigkeit, T die 
Umlaufzeit, t die laufende Zeit bedeuten. 

1. Fall: Konstante Winkelgeschwindigkeit. 

a) Bezogen auf den Brennpunkt: 

Mit 7 = 2 und der Integralformel?) 


2n 
27 


dp An 
J|il+kcosp yYi-r 
6 


. (5) 


erhält man nach (4a) und (1a) das bekannte Ergebnis?) 
2n 

u ep dp BERLIN ©) 

zul Prog Te 

0 


OR 


Der Mittelwert aller Entfernungen ist also gleich der 
kleinen Halbachse, wie auch in der Arbeit von Emers- 
leben gezeigt wird. 


b) Bezogen auf den Mittelpunkt: 
Mit »7 = 2 ergeben (4b) und (1b) 


2n 
u b dıy 
= m 
0% De 4 
0 
oder Ik 
& 2b 
ren BE). RN (8), 


wenn K. das vollständige elliptische Normalintegral 
1. Gattung bezeichnet. Unter Benutzung des Gauß- 


*) Siehe z.B. A. Sommerfeld: Atombau und Spektral- 
linien. Braunschweig: Vieweg 1924, 4. Aufl.,'S. 773. 

®) Siehe z.B. R. Rothe: Höhere Mathematik, Bd. II, 2. Aufl., 
5.80. Teubner 1931, 
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schen arithmetisch-geometrischen Mittels M gilt 
bekanntlich‘) ; 


>» moin u) 
[77 gt! a 


ergibt (8) als Mittelwert aller auf den Ellipsenmittel- 
punkt bezogenen Entfernungen bei konstanter Winkel- 
geschwindigkeit den einfachen Wert 


ab 


Ein Vergleich von (6) und (8) ergibt, da K> = 


A>rı I <k<I 


Da alle Mittelwerte in Abhängigkeit von der Exzen- 
trizität k bei konstanter großer Halbachse betrachtet 
werden sollen, schreiben wir für (6) 


für 


7, = all ka. (12) 
und für (8) 
ee 
&=— 1%? K(k) (13). 
Mit X =Y1-—k2 gilt®) 
RER LEN 20 (14), 
X-0] ER 


so daß 
END, 2a 4 
lim 9, = — limk’K(k) = — lim % In — =0 (15). 
ET 01 er) ( ) 0 k’ ( ) 
Für k=1 sind also beide Mittelwerte Null (s. Bild) 


2. Fall: Konstante Flächengeschwindigkeit. 


a) Bezogen auf den Brennpunkt. 
Diesen Fall behandelt Emersleben. Da bei kon- 
stanter Flächengeschwindigkeit 


ae Dit 
EN 
SPP=T 0...) 
gilt, ergeben (4a) und (la): 
5 an , 
RR dp am. 
a Zadar |(1-+ kcosp)® N 
0 
Für das letzte Integral findet Emersleben 
27 fen 
er. In RE (18 
ka ua 
{) 
so daß (17) nach einfacher Umformung in 
2 e 
a=al147): Du Dr 10) 


) 
übergeht. Dieser Mittelwert ist größer als a und 
nähert sich für k — 1 dem Grenzwert 3a/2, während 
die unter 1 betrachteten Mittelwerte sämtlich dem 
Wert Null zustreben. - 


b) Bezogen auf den Mittelpunkt. 
Wie unter a) erhält man 


2n 2n 
dy 


er U bh? \ 
tu NE ag 2 2 See ee 
ea ee he ie (20), 
0 N) EA 


*) Kommerell: Das Grenzgebiet der elementaren und höhe 
ren Mathematik. Leipzig: K.F. Köhler 1931. ’ “ 
5) Jahnke-Emde: Teubner 1938, 


Funktionen-Tafeln, 
3. Aufl., 8. 73, j Sr 


x 
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wofür auch wieder 
n,2 
dy 


DL SET ER EN 
z VaZ—Rrcosy j 


02 etz) 


geschrieben werden en Dieses Integral ist aber 
durch das elliptische Normalintegral 2. Gattung E 
ausdrückbar®), nämlich 

n/2 


dy E(k) 


me me en . 22 y 
„ Na—rroosyp IR ie 
so daß folgt 
ern Elk ERPT 
a N, bw 9=%.E) (2). 
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3. Fall: Konstante Bahngeschwindigkeit v. 
&) Bezogen auf den Brennpunkt. 


Nach (3a) ist 


T. Mi 2 
Ir 1 1 1 ee 
n=-7 ra: = fras _ allesamt (27). 
ö ö 0 


Da der Ellipsenumfang 
4aE(k)=vT 
ist, ergibt sich mit (1b) nach einfacher Umformung 


SAID NEREIETEN., 
AT IaE eg 


. (28) 


(1+ kcosp)°® 


BRUNS 


“u. SERIEN 


Man erhält also einen Wert, der ganz entsprechend 
gebaut ist, wie @, nach (8), Wie das Bild zeigt, fallen 
0, und ®, für kleine &2 praktisch zusammen. Diese 
bestätigen auch die Reiherientwicklungen: 


2 k® or 
nei 4-0.) 


3 .) | 


n ....(24), 
u= a(ı Be Fuss 


4 64 


die zeigen, daß d, für kleine k? sich nur um Glieder 
von der Größenordnung k? von 6, unterscheidet. 


Stets ist : 
IN E 
0,> 0, da grhı-e. a 12.628) 
und 
| im, = -ü (26) 
21 2 er EL ae EN TR} . 


9) Jahnke-Emde:l.c. 8.56. 


"Führt man in dem letzten Integral = g92ein und 


benutzt die Beziehung cose2%=1--2sin?#, so 
erhält man s 


V(1-+ k)?— 4ksin?# 
ker  llnear, rar ze 
Führt man die Abkürzungen ein 
4k A 
=k?< 
(FR ESS IR AN, (30) 
2k P D 
ar kılk=n<sI1l, 
so geht (29) über in 
n/2 A 
— kügi 
x a N Yl— k!sin? 9 49 (31). 


n=7m1M N (1— nsin2 9) 


u 


72 


TEEN, 


ee 


Re 


FR" 
£ 


Enge FR, a: 
I EN AR 
a ee 
. 6, es 


ca 
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En: 
Setzt man sind = u, so wird 
1 
1— kiu? 


T= ; 
HDy ..(32). 
x ? 


ao) — Hab) 


Wendet man auf die rationale Funktion unter dem 


Integralzeichen Partialbruchzerlegung an und führt 
die Abkürzung ein 


1 
du 


- ——d 
fü — nu2)’ /(1— 2) (1 — kiu?) 


v ...(33), 


so folgt 
(1—R)? 
1l+%k 


Die Integrale Z, lassen sich mit Hilfe einer Rekur- 
sionsformel?) auf Z,, Z,, Z-, zurückführen. Die Durch- 
führung ergibt: 


a a 
\  E(k) 


N 
(22,—(1—k)2,) . (34 


T,= 
a lb 1 le . (35). 
Nun ist aber 
du 
Gi k ..(36), ° 
e means" Sa en B 
Il—nu DER 
— K(k 
et 3 NER (37) 
+ 5 (Bl) — Rh) 
oder unter Benutzung von (30) “E 
1% 3% 
Z1= "SE Ki) + SEE) 2.88 
und 
N $ 
2: nn ‘2 ne. 
z Ju-au) Ya—u)(1—ru) (. (89. 
= II(n, k,) SE ' 


so daß (35) in 


das vollständige elliptische N ormalintegral, 3. Gattung, = 


2% 


rs - sla-nzm | 


er 


übergeht. Pe: t 
Die numerische Auswertung von (40) mit Hilfe der 
Legendreschen und Heumannschen?) Tafeln legte 
die Vermutung nahe, daß 7, den konstanten Wert «a 
hat. Der Beweis dieser Vermutung gelingt wie Be: 
Löst man die Legendresche Beziehung?) Ev; 
z) m 
2 E(k)—(1+k) E(k,) 2 

1-2. ve? 


’) Serret-Scheffers Lehrbuch der Ditferenti nd 
tegralrechnung, Bd. II, 4. u. 5. Aufl., S. 93. Teubner 191. N: 

Carl Heumann: Tables of complet elliptie. äkpkeaIn. 
J. Mathematics Physics. Vol.XX, Nr.2 19a). A. 

?) Enneper: Elliptische Funktionen. Fe 1876, 
Formel (21). 


K(k) = 


nach B(k,) auf und führt es in (40) ein, so folgt 


1%: 
173, = a+ +) E(k) ae k)II(n, k,) . (42). 
— (1+k)K(k)} 
Da nach der Gaußschen Transformation 
5 et 3 
rn IE BikY a ee (43), 
hat man 
SE 
z,=0+— 5 a - {(1L— k) II (n, k,) — K(k,)} (4). 
Ferner ergibt (30) 22 
ö 2Kk, Bu 
I-k-gug n=l-h 2.) 
mit ; 
k; Zan vı 7 k? ’ 
so daß (44) in 
1% > 
ln Fpu-H, &)— Kiky)| (40) 


„nämlich 


so folgt nach Rinfacher.  Umformung i 


mit 


ar auch 
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übergeht. Benutzt man schließlich eine ebenfalls aut 
Legendre zurückgehende Beziehung!'), die es ge- 
stattet, den Parameter eines elliptischen Integrals 
3. Gattung auf einen anderen zu transformiren, 


— K(k) + 08? & IT (k? sin? a, k) \ - 


k®sinta , ee BER 
1— k?sin?a 1 — k?sin?a’ )=- 
und setzt in ihr tg’«—= 1%’, so wird 
 k’?sin?a k’ 
’ 2 ZZ m 0 
ET I k2sn2a. , 1EPRT (48) 
-  k2cos?& 
EN WE EL 
| Kein?a=  — ya sing 1—W...(4) 
und (47) ergibt Dr 
— Lk’ hy 
iR iR a k ‚k)= ” (50). 


Denkt man sich für k k geschrieben, so NE (0) 
das Verschwinden der eckige n Klammer in (46) any 
womit r man den Beberehett einfachen Wert se 

eh=e Be ee 2 . 
- erhält. 


. .b) Bezogen auf den re re 
ae in (27) gilt net DR = & KEN Fee 


A ferer er dy. RER er = 


ale man a ein führt die Ditterentiation aus OH 


FRE By, BR 
3 ur Baal a “ 


a 


1 De 
« 
c 3 * \ z k & 
N 2% ’2 
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Führt man dieses elliptische Integral auf die Normal- 


form zurück unter Benutzung des von (32) nach (40) 

führenden Er so erhält man 
re r 
= er: oe)... (m 


oder . 


. (58) 
mit 
M=R2—=%2), KB=LI—R:. 

Der Grenzwert dieses Ausdruckes ‚für k—1 ergibt 

sich zu a/2, da wegen (14) 
lim Back) ten (59) 
.e 
und 
lim B(k) = lim BE) Er nr k00): 
k>1 k>1 


Die numerisch festgestellte Größenbeziehung 0,<z 
<’0, wird auch für hinreichend kleine k? durch die 
Reihenentwicklung 

= Eee 
ao m- 


im Vergleich zu (24) bestätigt. 
Alle betrachteten Mittelwerte sind auf dem Bild 
für «= 1 in Abhängigkeit von % dargestellt. 


Braunschweig, A. Fricke. 


) k?<1 . (61) 


Fehlerausgleichung bei Funktionen in Para- 
meterdarstellung. 


Während das Verfahren zur Fehlerausgleichung 
explizit gegebener Funktionen allgemein bekannt ist, 
findet man dieses Problem für Funktionen in Para- 
meterform nicht dargestellt. Als Verallgemeinerung 
der bisher üblichen Aufgabe soll es im folgenden be- 
handelt werden. Die Anregung dazu ergab sich aus 
Untersuchungen über das Dampfgleichgewicht ver- 
schiedenatomigen Schwefels im hiesigen Phys.-Chem. 
Institut. 

Der Zusammenhang zweier Variablen x und y sei 
in Parameterdarstellung gegeben 


fl,  y=olbW).... 
mit dem laufenden Parameter it und der Konstanten u, 
deren, Wert aus zwei Er Meßreihen. 


BERURX N. „An Mi 
Erz. RER) 


zu bestimmen ist. 
Zunächst möge zur Kennzeichnung des Verfahrens 


der bekannte Sachverhalt skizziert werden für den 


Fall, daß der Parameter it eliminiert werden kann. 
Man erhält dann den Zusammenhang in der Form 


\ UA Rs eis ers ern (3). 


Die Konstante u bestimmt man nach dem Prinzip 


des kleinsten Fehlerquadrats, indem man die Ordi- 


Er natenabweichung von Meßwert und gesuchter Kurve 
bildet und ihre Quadratsumme zum Minimum macht: 


nn Sn, uf = Min RR 
ö ENTER i 
Die Meßwerte X, werden dabei als fehlerfrei an- 
;ehen. 
' Dieses im allgemeinen nichtlineare Problem lineari- 
"siert, man mit Hilfe eines geeigneten Näherungswertes x 
HR ‚erhält mit u =ü + Öu bei Beschränkung auf das 
are Glied der Korrektur öu eine einfache Ersatz- 
= Ex En: * 


es A Er 2 re a 
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Mit den Abkürzungen 


Yv (Xr, u), Zu u (X, u), 
indem man die verschiedenen Abszissen durch den 
Index v» kennzeichnet, ergibt sich 


hiX,u)=w-+ U gu 


und damit an Stelle von (4) % Ba; Form 


£ oyv z Y = 
DE usa nin du = Min, , » , (5), 


deren Normalgleichung linear in der Korrektur du 
wird. Das Verfahren wird sodann mit dem ver- 
besserten Näherungswerte @ = u + öu wiederholt und 
solange fortgesetzt, bis die verlangte Genauigkeit er- 
reicht ist. 

Wenn nun die Elimination des Parameters t in (1) 
nicht möglich ist, so liegt die Aufgabe vor, die 
Quadratsumme der Ordinatendifferenzen’ zum Mini- 
mum zu machen 


[Hr —gli,u)®=Min .... (6 
v=1 


unter der Nebenbedingung 


ED ER NRZ); 


die zum Ausdruck bringt, daß diejenigen Ordinaten 
zu wählen sind, die zu den festen «-Werten der Meß- 
punkte gehören. Infolge der Nebenbedingung hängen 
die Parameterwerte {, in komplizierter Weise von der 
gesuchten Größe u ab, so daß man im allgemeinen 
linearisieren, muß. 

Wir denken uns mit einem geeigneten % eine Nähe- 
rungskurve ermittelt, die infolge der Parameter- 
darstellung nach t beziffert ist, und darauf diejenigen 


Punkte &, y, mit ihren, Bezifferungen {, markiert, die 
zu den Meßpunkten X, Y, gehören, und zwar so, daß 
I = Ay wird, 


ACER RN (8). 
Mit diesen Näherungswerten % und 1%, schreiben wir 
u=u+öu, hb=b+öb 


und entwickeln die in (6) und (7) auftretenden Funk- 
tionen bis zu den linearen Gliedern 


Bee Chiya Of 

flv, u) = fl, u) + ER a)öl, + ER a)öu (9), 
® du - i\ 

g9(t, u) = g(tv, u) + Eli u) 6 + IA a)öu (10). 


Wegen (7) und (8) entfallen die ersten beiden Glieder 
in (9), so daß sich aus dieser Gleichung öt, durch öu 
ausdrücken läßt. In (10) eingesetzt ergibt sich 


Oyv , 02) ou. 


ou "ou 
wobei zur Abkürzung gesetzt ist: 


9 (lv, u) = yr+[ 


er ER oyv _0g - 
RR g(to; u) E ou zn (h, u) 
097% 
BLAlre 
BE RERSEN 90, Of - = 
v Br FAZ u). 


of - _ 
3. u) ’ 
Damit ist in der linearisierten Form der Parameter 
eliminiert und es ergibt sich an Stelle von (6), (7) eine 
6) RT Formulierung 
Oy» ‚0x 
I [y,_ Unc- [er — Yr du ee) gu| = Min. (I 


=. l 


186 


Kleine Mitteilungen 


Z. angew. Math. Mech. 3 


die in der üblichen Weise weiterbehandelt werden 
kann. Weitere Konstante neben u lassen sich analog 
berücksichtigen. 

Ähnliche Darstellungen ergeben sich, wenn man 
nicht wie hier die.Ordinatenabstände, sondern die 
Abszissen- oder die Normalenabstände oder auch die 
{-Differenzen in (1) mit (2) unmittelbar ausgleicht. 


Hannover. K. Jaeckel. 


Vektoren auf der Kugelfläche. 

Zu den Vektoren der Ebene kann man ein inter- 
essantes Analogon von gekrümmten Vektoren auf der 
Kugelfläche angeben, Verbinden wir nämlich mit dem 
Bogen eines Großkreises der Kugel gleichzeitig eine 
Richtungsangabe, indem wir den Betrag des Bogens 
durch c, seinen Durchlaufungssinn mit dem Rich- 
tungsfaktor c® bezeichnen, so gilt für den ge- 
richteten Bogen c, den wirals Kugelvektor 
bezeichnen wollen, 

ae (). 


In bezug auf den Bogen von der Größe c spielt also der 
Richtungsfaktor c® lediglich die Rolle des positiven 
oder negativen Vorzeichens. Analog den goniometri- 
schen Definitionen 


Gcile 


sin (—2) = — sin?2; c08 (—2) = 008% 


definieren wir für Kugelvektoren 
(2). 


Wir definieren ferner in der Art gewöhnlicher Vek- 
toren für ein sphärisches Polygon die (nichtkommu- 
tative) Komponentenzerlegung (Bild 1) 


Behrens 


(RE 2 
\J 


sine = sin cc! = (Psin c; 
COS C= 6086! = 008 6 


PR T ES 


Bild 1. 


und verlangen, daß auf Summen von Kugelvektoren 
die Additionstheoreme der elementaren Goniometrie 
angewandt werden können. Hierbei treten Produkte 
der oben eingeführten Riehtungsfaktoren auf, die wir 
Richtungsprodukte nennen wollen. Die Be- 
deutung solcher Richtungsprodukte von der Form 
” (a°. 6°) ergibt sich aus einfachen Sonderfällen. 


1 
a) Aus = cH36 d.h. 


08 c= 008 - cos 
2 2 
oder 
c \ {H 
coSC—= 008? 55 — sin? (cc) 


folgt wegen 


sofort 
(lose Zn were). 


b) In einem gleichschenkligen sphärischen Dreieck 
(Bild 2) mit’ den Seiten «= b=. 90° gilt 
c= 90° a + 906%, d.h. 
co8 c = 0082 90° — sin? 90% (a0 59) , | 
Folglich muß sein 
(a 6%) = — 008c, 
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d.h. gleich dem cos des von den positiven Richtungen 
von’a® und b° eingeschlossenen Winkels (180° — c), hat 
also dieselbe Bedeutung wie in (4). 

Wir definieren also allgemein: Das Richtungs- 
produkt (a°-bP) ist gleich dem cos des von den 
positiven Richtungen eingeschlossenen Winkels. Es 
ist kommutativ: 


(a° 50) = (b9 - 0%) = cos (ad, 6), . . . ). 


( C $ 
P14 BR 


a 


Bild 2. 


Für ein beliebiges Dreieck aus Kugelvektoren (Bild1) 
ist c=a-+b oder 
c08c= 0084 c0osb — sin asin b (a0 6°), 
008 0= cosacosb — sina sinb cosYy. 


Das ist der bekannte Seiten-cos-Satz der sphärischen 
Trigonometrie. 


d.h, 


Das soeben betrachtete Richtungsprodukt ergab, 
sich unmittelbar als Folge der trigonometrischen cos- 
Operation. Wir können eine ähnliche Festsetzung 
auch für die sin-Operation treffen: 


a) Aus c= = c+ = c ergibt sich 


esinc=2 sin >. 0085 
oder nach Richtungsmultiplikation mit c’ wegen 
(ce? 0) = 1 richtig 


sinc= 2sin x cos 2 
— In ee 
2 2 


b) Im Dreieck von Bild 2 ist 
c? sin c= a? sin 90° cos 90° + 6° sin 90° cos 90°. 

Da die rechte Seite 0 ist, bleibt diese Beziehung nur 
richtig, wenn die sin-Operation noch mit einer sofort 
anschließenden Richtungsmultiplikation mit a® oder 
b° verbunden wird, da ja (c?-a9) = (c°-b%)= 0. Da 
aber ferner (c? - c0) = list, darf nach der sin-Operation 
nicht mit c° riehtungsmultipliziert werden. 

Wir gewinnen somit den Satz: 

Nach der sin-Operation muß eine Richtungsmulti- 
plikation mit einem der Richtungsfaktoren der' Kom- 
ponentenvektoren angeschlossen werden, wenn nicht 
eine andere, weiter unten zu besprechende Drehmulti- 
plikation vorangeht. 

In einem beliebigen Vektordreieck (Bild 1) mit 
ce= a + b ist zunächst 
 esine= a'sinacosb +bPcosa sind. 
Eine Riehtungsmultiplikation mit a® bzw. b° liefert 
wegen 
(ad cd)= cosß; (aP-bP) = cosyY; 
den bekannten sin-cos-Satz: 
sin c cos ß = sin a c08b + cosa sin b cos y; 
sin c 608% = cosa sinb-+ sin«a cosb cosy. 


(69 00) = cosa 


Die weiteren durch zyklische Vertauschung hierausent- 


stehenden Formeln ergeben sich auch in analoger Weise 
aus den Darstellungen er 


=e-bi  der-a, 


Um auch die Winkel zwischen den Richtungen zweier 
Kugelvektoren erfassen zu können, führen wir ein wei- 


teres Produkt ein, welches wir aus sogleich zu erör- 


ternden Gründen Drehprodukt nennen wollen, ) 3 & 


“ J we 
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Drehen wir den Kugelvektor a (Bild 3) im positiven 
Winkelsinne in die Richtung von b, so beschreibt der 
| positive Pol (d.i. der Pol, von dem aus die Rishtung 
von a im positiven Winkelsinn erscheint) den Groß- 
kreisbogen €, der somit im positiven Pol von b endet. 
Es ist also 


C=YEN.deh, 
Wir defioieren zweckmäßig das Drehprodukt 
a NT N EEE (7), 


damit bei senkrechten Vektoren a und b das Drehpro- 
dukt ihrer Richtungen die Richtung des Polarvek- 
tors € ergibt, während das Drehprodukt gleicher 
Richtungen verschwindet: 


{a Se, — 102.1 22a): 
Das Drehprodukt ist alternativ: 
(ax) = = (b'Xa9)'. 2%. (76). 


Zu einem sphärischen Dreieck (Bildl) a+b=c exi- 
stiert ersichtlich ein Polardreieck 1 +b= derart, 
daß 
(ax b) =R'siny; 
(a°x c%) = 6° sin ß; 
L- (0x6) = Asina; (EC x6) = asina; 
(a- 6%) = cos c; (6°.70) = cosa; 
(€. a°)= cos b. 
Dies ist der bekannte Ausdruck dafür, daß im Dreieck 
und Polardreieck Seiten und Winkel vertauscht sind. 


(6°x 60) = (sin c; 
(ax ce) = b° sin d; 


(8). 


Bild 3, 


Für ein gegebenes Dreieck c=a + b auf der Kugel 
wenden wir die sin-Operation an und erhalten 


ce sinc= a’sinacosb + b° cosa sin b 

oder nach Drehmultiplikation mit a® oder b° 
(e?x a®) sin c= (ax a°) sin a cosb+(b°x a°) cosa sinb; 
bzw. hi 

(x b°) sin c=(a°x b°)sin a cos b+-(b°x b°) cos a sin b, 
d.h. nach (8) und (7a, c) 

— b° sin csinß= — (° cosasin bsin y; 

bzw. a’'sinesing= (sinacosbsiny. 


Eine Richtungsmultiplikation mit € führt diese beiden 
Gleichungen mit Rücksicht auf (8) über in 


" sincsinß=sinbsiny, bzw. sincsina = sinasin y 


Das ist der bekannte sin-Satz. — Man überzeugt sich 
Ar auch hier sehr leicht davon, daß ebenfalls die Dreh- 
-  multiplikation stets nur von der Kompo- 
nentenseite her erfolgen darf, so daß der 
oben ausgesprochene Satz folgende neue Fassung er- 

‚hält: ; 


plikation oder eine Drehmultiplikation mit nachfol- 
. gender Richtungsmultiplikation angeschlossen werden, 
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wobei alle Multiplikationen nur von der Komponen- 
tenseits her erfolgen dürfen. 
Wendet man die bisherigen Entwieklungen auf das 


Polardreieck a+b= c an, so folgen die bekannten 
Polarformeln, auf deren Ableitung hier verzichtet 
werden kann, da sie dem Vorstehenden analog erfolgt. 

Die betrachteten Zusammenhänge lassen sich mit 
bestimmten Einschränkungen auch auf sphärische 
Polygone ausdehnen. Wir behandeln als Beispiel ein 
sphärisches Viereck (Bild4), bei dem die Seiten a, 
b, c und die Winkel« und ß gegeben sind. Wir schreiben 


dD=a+b-+e. 


Hier ist wegen der nicht-kommutativen Eigenschaft 
der Summe die richtige Reihenfolge der Vektoren wich- 
tig. Aus den beiden Anordnungen 


dD=(a+b)+c 
b=a+(+e) 


bzw. 


ergibt sich sofort 


cos d= cos(a + b)cose — sin(a-+b)sinc 
bzw. 
cos d= 0osa cos(b + c) — sinasin (b + c). 
In beiden Fällen liefert die weitere Zerlegung wegen 
(a -69%) = c0sa; (69-9) = cos: 
c08 d= 608 4 608 b 60Ss c— sin a sin b COS € COS& 
— sin a cos b sin c (a + c®) 
— cosa sin bsin c cos ß 


(9). 


> 


Bild 4, 


Hier besteht noch die Frage nach der Bedeutung des 
Richtungsprodukts (a®-c°%) der nicht zusammen- 
stoßenden Vektoren a und c. 

Eine Hilfsbetrachtung im Viereck AABB, in dem 
AA und BB gleich 90° gemacht worden sind, liefert 
für AB = b sofort 

| b= 90° 0° +65 4900 c0, 


dh, 
cosb= — (ac) cosb 
oder 
Be cos b 
(a? » c0) RR (10) 


Wird 5 = 0, a. h. geht das Viereck in ein Dreieck über, 
so wird d= 1800 — Xl(a,c) und cosb=]; es ent- 
steht also unsere alte Definition (5). Man erkennt aus 
der Figur nun leicht folgendes: Errichtet man in A 
und B die zu den Polen führenden sphärischen Senk- 
rechten auf a und c, so schneiden sich diese unter einem 
Winkel o, Die Drehung des Dreiecks AM Ä um den 
Schnittpunkt M führt den Punkt Ä nach einer Drehung 
um den Winkel bin 3 über; dann fallen aber die Rich- 
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tungen von AM und BM zusammen. Also ist 180° — o 
— b, so daß (10) übergeht in 


(MER 


.(ul). 


hr, Vale Kot 


cos b 


Durch eine vektorische Rechnung in der obigen Art 
oder einfacher durch Anwendung des sphärischen Win- 
kel-cos-Satzes, der ja oben vextorisch herleitbar war, 
findet man im Dreieck ABM 
co80o=—sinasinß + cosacosßcosb. 
Führt man dies in Verbindung mit (11) in (10) ein, 
so ergibt sich der endgültige Ausdruck für cos.d, den 
wir auch nach längerer umständlicher Rechnung durch 
Benutzung der bekannten elementaren sphärischen 
Sätze erhalten hätten. — Um den Winkel ö zu erhalten, 
wenden wir die sin-Operation aut d=a-+(b-+c)an 
und erhalten nach Richtungsmultiplikation mit c® 
sin d cos ö = sin a cos b cos € (a® - c®) y 
+ cos a sin b cos c c0sß + cos a cosbsin c 
— sin asin bsin 6 c0sa, 


BUCHBESPRECHUNGEN 


Pickert, Einführungin diehöhere Al- 
gebra. (Studia Mathematica — Math. Lehrbücher 
Bd. VII.) 298 S. Göttingen 1951. Vandenhoeck & 
Ruprecht. Preis brosch. 12,80; geb. 14,80 DM. 

Eingeordnet in den Begriff der algebraischen Struk- 
tur, werden die wichtigsten Grundbegriffe der moder- 
nen Algebra entwickelt und bis zu den Bewertungen 
von Körpern fortgeführt. Geschickt ist die Zusammen- 
stellung der einem jeden Paragraphen hinzugefügten 
Aufgaben, die nicht nur den Zweck der Einübung 
des jeweils voraufgegangenen Stoffes dienen, sondern. 
dem Laser auch zusätzliches Wissen und, die Einsicht 
in weitere Zusammenhänge vermitteln sollen. So ent- 
steht im ganzen ein klares und trotz des knappen 
Raumes abgerundetes Bild von der modernen alge- 
braischen Forschung. 

Nur in einem erscheint mir die Auffassung des Ver- 
fassers reichlich optimistisch: daß der Studierende 
nach dem zweiten Semester bereits der geeignete und 
berufene Leser dieses Werkes sei. Für ihn besteht die 
Gefahr, den hohen Allgemeinheitsgrad als bloßen For- 
malismus zu nehmen, weil er nicht zu erkennen ver- 
mag, daß es sich um die Zwangsläufigkeit geschichtlich 
notwendig fortschreitenden Formalisierens handelt. 
Erst auf der Grundlage bestimmter algebraischer 
Kenntnisse kann der Wert des vorliegenden Buches 
voll ausgeschöpft werden. ö 


- Dresden. Draeger. 


Mitteilungen aus dem Max-Planck-Institut für 
Strömungsforschung. Unter Mitwirkung von L. 
Prandtl, W. Tollmien, G. Vogelpohl. 
Herausgegeben von A. Betz. Selbstverlag Max- 
Planck-Institut für Strömungsforschung Göttingen. 
Nr.1. A. Betz, Zur Einführung. — J.Rotta, 
Über die Theorie der turbulenten Grenzschichten. 
548. mit 10Abb. Göttingen 1950. Preis brosch. 
4,— DM. 

Nr.2. H. Himmelskamp, Profilunter- 
suchungen an einem umlaufenden Propeller. 38 S. 
mit 59 Abb. Göttingen 1950. Preis brosch. 8,— DM. 

Um eine Möglichkeit zu schaffen, geeignete Arbeiten 
größeren Umfanges in der für spezielle Interessenten 
erwünschten Ausführlichkeit zu veröffentlichen, ohne 
die bestehenden Zeitschriften allzusehr zu überlasten, 
hat sich der Direktor des Max-Planck-Instituts für 
Strömungsforschung, Prof. Betz, entschlossen, 
solche neuere Arbeiten aus seinem Institut, sowie 
einzelne ältere, nicht zugängliche Arbeiten in dem 
preiswerten Rotaprintverfahren zu drucken und m 
zwangloser Folge herauszugeben. 


Bd.31 Nr.6 Juui1l 090 
was mit Rücksieht auf (11) nach längerer Rechnung 
auch direkt hergeleitet werden könnte. Ebenso er- 
gibt sich ein Ausdruck für sin d cosy durch Anwendung 

der sin-Operation auf d= (a-+b)+c und nachfol- 
gender Richtungsmultiplikation mit a®. Beider Berech- 
nung von ö und y ist also zu beachten, daß die Klam- 
mern nd=a-+b--c so gesetzt werden, daß die 
dem gesuchten Winkel nicht,anlie- 
genden Seiten in Klammern ein- 
geschlossen werden. — Ähnlich läßt sich ver- 
fahren bei beliebigen Polygonen. 


Die Verwendung der beschriebenen Kugelvektoren 
ist auf Fälle der vorstehenden Art keineswegs be- 
schränkt. Sie erweisen ihre Nützlichkeit z.B. auch 
bei Fragen der analytischen sphärischen Geometrie. 
Auf diese Fragen wird in einer späteren Arbeit geson- 
dert eingegangen werden. 


Elmshorn, H. Athen. 


Die Arbeit von Rotta gibt an Hand früherer und 
eigener neuer Überlegungen und Rechnungen, sowie 
vorhandener Messungen einen schönen Überblick über 
das, was man heute über die turbulenten Grenz- 
schichten aussagen kann, und welche Meßreihen er- 
forderlich erscheinen, um Näherungsverfahren für die 
Berechnung turbulenter Grenzschichten zu entwickeln. 

Im besonderen wird nachgewiesen, daß auch bei 
turbulenten Grenzschichten, wie bei den laminaren, 
„ähnliche‘‘ Lösungen existieren, bei denen das Ge- 
schwindigkeitsprofil längs der Wand nur affin ver- 
zerrt wird. Für diese ähnlichen Lösungen werden 
dann unter Zugrundelegung vorliegenden Versuchs- 
materials Näherungsrechnungen durchgeführt. 

In der Praxis hat sich gezeigt, daß eine Übertragung 
von Eigenschaften feststehender Gitter auf Laufrad- 
gitter recht unsichere Ergebnisse liefert. Insbesondere 
treten am umlaufenden Profil erhöhte Auftriebsbei- 
werte auf, die durch die Rotation bedingt sein müssen. 
Um hier zuverlässige Rechnungsgrundlagen zu 
schaffen, hat Himmelskamp Druckverteilungs- 
messungen an einem Blatt eines vereinfachten Ge- 
bläselaufrades angestellt und sie mit entsprechenden 
Ergebnissen in ebener Strömung verglichen. Außer- 
dem wurde an einigen Punkten des Blattes die Rei- 
bungsschicht ausgemessen und damit ein näherer Ein- 
blick in die Vorgänge einer rotierenden Reibungs- 
schicht gewonnen. 

Diese, in zahlreichen Kurventafeln niedergelegten 
Messungsergebnisse dürfen als ein wertvoller experi- 
menteller Beitrag zur Auffüllung der Lücke in den 
Rechnungsgrundlagen bezeichnet werden, die noch 
für die Berechnung von Laufradgittern besteht. 


Weilburga.d. Lahn. * L. Schiller. Si 


Dr. phil, Dr. Ing. e.h. Werner Schmeidler (vorm. 
o. Prof. a. d. Technischen Hochschule Berlin), 
Integralgleichungen mit Anwen-. 
dungenin Physik und Technik. Bd.I: 
Lineare Integralgleichungen. (Mathematik und ihre 
Anwendungen in Physik und Technik, ReiheA, 
Bd. 22.) VII u. 6118. Leipzig 1950. Akademische 
Verlagsgesellschaft. Preis geb. 3,— DM. 

Das vorliegende Werk ist die erste deutsche lehr 
buchmäßige Gesamtdarstellung der Theorie de 
Integralgleichungen. Es füllt in mustergültiger Weise 
die bisher zwischen den bekannten guten Ein 
rungen und der Originalliteratur klaffende Lüc 

Charakteristisch für das Buch ist die durchla 
‘Benutzung der Hilbertschen Methode ( 
ordinatenfunktionen — selbstverstä ıdlichi 

N Kay ® 
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Fassung. Dadurch wird der umfangreiche Stoff ein- 
heitlich zusammengefaßt und gleichzeitig eine orga- 
nische Verbindung zwischen den Integralgleichungen 
1., 2. und 3. Art hergestellt. Der Verfasser konnte 
hierbei — insbesondere für die Integralgleichungen 
1. Art — auch neue eigene Ergebnisse benutzen. In 
und neben dieser grundlegenden Darstellung der 
Theorie finden auch die wichtigsten andersartigen 
Ansätze und Verfahren ergänzend und abrundend 
ihren Platz. Einen Einblick in den reichhaltigen 
Inhalt sollen einige Stichworte aus dem Inhalts- 
verzeichnis vermitteln. 

Theorie der Integralgleichungen 1. Art, dazu spe- 
zielle Formen. wie Fouriertransformationen und Ver- 
wandtes, Volterrasche Gleichungen, Beispiele 
und praktische Lösungsverfahren. Theorie der Inte- 
gralgleichungen 2. Art, Volterrascher Typ, qua- 
dratisch integrable Kerne einschl. der Fredholm- 
schen Theorie, die Extremwerteigenschaften von 
Eigenwerten, der Zusammenhang mit Differential- 
gleichungen, Lösungsverfahren, singuläre Kerne 
einschl. der Carlemanschen Theorie, Integral- 
gleichungen 3. Art. Ein Anhang behandelt unend- 
liche Gleichungssysteme. 

Etwa 100 sorgfältig durchgeführte Beispiele be- 
leben nicht nur die theoretischen Untersuchungen, 
sondern werden zusammen mit zahlreichen Abschnit- 
ten über numerische Lösungsmethoden auch den an- 
gewandten Mathematiker bei seiner Arbeit unter- 
stützen. 

In der Darstellung bestehen zahlreiche Parallelen 
zu den nun schon bekannten ‚‚Vorträgen über Deter- 
minanten und Matrizen in Physik und Technik‘‘ des 
Verfassers [Besprechung in dieser Zr Bd. 29 (1949) 
S. 287]. Ebenso wie dort werden nur die bei Hoch- 
schulingenieuren üblichen Vorkenntnisse und einige 
Übung im mathematischen Denken vorausgesetzt. 
Notwendige Ergänzungen hierzu (Lebesguesche 
und Stieltjessche Integrale, über Orthogonal- 
systeme und Bilinearformen) bietet ein besonderer 
Anhang. 

Mit seinem reichhaltigen, Inhalt wird das Buch in 
. kürzester Zeit jedem, der sich mit Integralgleichungen 
befassen. will, unentbehrlich sein. Seine Benutzung 
als Nachschlagewerk wird durch die übersichtliche 
Zusammenfassung der Einzelresultate zu Sätzen, 
durch ein Verzeichnis aller behandelten Beispiele so- 
wie durch einen bis etwa 1944 reichenden Literatur- 
nachweis erleichtert. 


Dresden. N.J.Lehmann. 


Dr. J. Horn (ehem. o. Prof. a. d. Techn. Hochschule 
Darmstadt), Partielle Differentialgleichungen. 
(Göschens Lehrbücherei, Gr. I, Bd. 14.) Vierte unver- 
änderte Auflage, 228 S. m. 8 Abb. Berlin 1949. Walter 
de Gruyter & Co. Preis 14,— DM. re 


. Allgemein wird begrüßt werden, daß das seit längerer 
Zeit vergriffene Buch des vor einigen, Jahren ver- 
_storbenen Prof. Horn wieder zu haben ist. Die Neu- 
auflage ist ein unveränderter Abdruck der dritten, die 
wiederum bis auf das ergänzte Literaturverzeichnis im 
Wesentlichen mit der 1929 erschienenen zweiten. über- 
einstimmte. Diese wurde in’ Bd. 10, $. 421/22 dieser 
Zeitschrift von Iglisch besprochen. Das didaktisch 
geschickt aufgebaute Buch wird wie bisher gern benutzt 
werden. 


“Dresden. Willers. 


E 


5 Die vor 
R Een ;t wieder einer gründlichen Überarbeitung 
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unterzogen; manches Neue ist hinzugefügt und, um 
den Umfang nicht wesentlich zu ändern, anderes 
fortgelassen. Die große Fülle des außerordentlich an- 
regenden Stoffes ordnet es in drei Abschnitte ‚‚Allerlei 
Unterhaltungsmathematik‘‘, ‚Von, den Zahlen‘ und 
‚‚Von den geometrischen Formen‘, und bringt darin 
das Interessanteste aus einer Reihe von Einzel- 
veröffentlichungen aber auch viel Neues. Das Buch 
wird auch in der neuen Auflage viele Freunde finden; 
insbesondere werden die Lehrer sein Wiedererscheinen 
begrüßen, da sie darin Vieles finden, das ihnen die 
Möglichkeit gibt, ihren Unterricht zu beleben und 
interessant zu gestalten. 


Dresden. Willers. 

Josef Lense (o. Prof. a. d. Techn. Hochschule 
München, Vom Wesen der Mathematik 
und ihren Grundlagen. 688. München 
1949. Leibniz-Verlag. 

Eine faßlich und klar geschriebene Broschüre, die 
auf geringem Raume einen guten Einblick vermittelt. 
Das gilt in besonderem Maße von dem kurzen Ab- 
schnitt über den mathematischen Intuitionismus, 
dessen Wesen sehr anschaulich klar gemacht wird. 

Am ausführlichsten werden Fragen der Mengenlehre 
und der Logistik behandelt. Das ist deshalb bedauer- 
lich, weildadurch aktuellere das Wesen und die Grund- 
lagen der Mathematik betr. Probleme überhaupt nicht 
zur Sprache kommen. Der Abschnitt über das An- 
wendbarkeitsproblem ist zu inhaltsarm, und die 
Wahrscheinlichkeit, die dabei gestreift wird, ist gerade 
heute auch ein intern mathematisches Problem. 

Dresden. Draeger. 


Dr. H. Küstner, Fünfstellige Logarithmen der 
natürlichen Zahlen und der Winkelfunktionen bei dezi- 
mal geteiltem Altgrad. VI-+157 S. Berlin/Leipzig 
1950. Verlag Volk und Wissen. 

Die vorliegenden fünfstelligen Tafeln, sind für den 
Gebrauch an Ingenieurschulen bestimmt. Sie enthalten 
daher außer den üblichen übersichtlich angeordneten 
mathematischen Tafeln umfangreiche physikalische und 
vor allen technische Tabellen für Maschinenbau, Hoch- 
und Tiefbau und Elektrotechnik. Wünschenswert wäre, 
daß in einer Neuauflage die Bereiche, in denen lineare 
Interpolation nicht erlaubt ist, genau bezeichnet 
würden. 


Dresden. Willers. 


Dr.-Ing., Dr.-Ing. e. h. R. Grammel (Prof. a. d. 
Techn. Hochschule Stuttgart). Der Kreisel, seine 
Theorie und Anwendungen. Zweite, neu- 
bearbeitete Auflage. Zweiter Band. Die Anwen- 
dungen des Kreisels. VI + 268 S. mit 133 Abb. 
Berlin-Göttingen-Heidelberg 1950. Springer-Verlag. 
Preis: geh. 30,— DM., geb. 33,— DM. 

Es ist außerordentlich erfreulich, daß so bald nach 
Erscheinen des ersten Bandes der zweite das Werk 
zum Abschluß bringt. Aufgabe dieses Bandes ist, 
zu zeigen, wie die im ersten abgeleiteten Sätze bei 
Behandlung des Kreisels in Maschinen bzw. des 
Kreisels als Maschine angewendet werden. Dieser 
Band ist vielleicht noch stärker umgearbeitet als der 
erste. Er berücksichtigt die Entwicklung der Kreisel- 
geräte bis zur Gegenwart, und wenn auch nicht alle 
Konstruktionen im einzelnen darin beschrieben wer- 
den — das kann nicht die Aufgabe eines solchen 
Buches sein —, so hat doch der Verf. alle wesentlichen 
Anwendungsgebiete berücksichtigt. Er teilt den Stoff 
in drei Abschnitte. Der erste behandelt die Kreisel- 
wirkung der Radsätze z. B. in Kollergängen und 
Pendelmühlen, bei stationär laufenden Radsätzen 
— insbesondere werden hier die kritischen Drehzahlen 
von. Rotoren untersucht — und in Fahrzeugen, 
Schiffen und Flugzeugen. Der zweite, umfangreichste 
Abschnitt ist dem Hauptanwendungsgebiet des Krei- 
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sels, den Kreiselgeräten, gewidmet, bei denen in der 
Hauptsache der schnellaufende symmetrische Kreisel 
benutzt wird. Zunächst werden Geräte mit Kompaß- 
kreisel und besonders der Kompaßkreisel selbst be- 
trachtet, ferner der als künstlicher Horizont benutzte 
Pendelkreisel mit vertikaler Achse und der als Kurs- 
anzeiger dienende Wende- und Lagekreisel mit mehr 
oder weniger waagerechter Achse und zum Abschluß 
Kreisel, die als Differentiatoren, Integratoren oder 
Regler dienen. In diesem Abschnitt wird auch be- 
sonders auf die Störungsfehler eingegangen. Der 
letzte kürzere Abschnitt beschäftigt sich mit unmittel- 
baren Stabilisatoren, also mit Richtkreiseln, wie es 
Planeten und rotierende Langgeschosse sind, mit 
Stützkreisen, durch deren Einbau’ labile Körper 
stabil werden (Einschienenbahn) und mit Dämpfungs- 
kreiseln, die die Schwingungen an sich stabiler 
Systeme abschwächen (Schiffskreisel). | 


Die beiden Bände des Werkes geben eine gut lesbare, 
umfassende Darstellung, die in gleicher Weise Theorie 
und Anwendung eines der interessantesten Gebiete 
der Mechanik berücksichtigt. Möge das Buch, das 
sich in gleicher Weise an Mathematiker, Physiker 
und Ingenieure wendet, von vielen durchgearbeitet 
werden. 


Dresden. 


Willers. 


Tables d’Interöts et d’Annuit6s editees par 
leCredit communal de Belgique. 163 S. Brüssel 
1950. 


Aus Anlaß seines neunzigjährigen Bestehens hat 
der Credit communal de Belgique diese Tafeln heraus- 
gegeben, die leider im Handel nicht zu haben sind. 
Die zum Teil auf 19 Stellen mit Maschinen berechneten 
Werte sind auf acht Stellen abgerundet hier abge- 
druckt und zwar für einen Zinsfuß von p—= 2 bis 8% 
mit einem Argumentschritt von 0,05%, so daß sich 
allgemein eine Interpolation erübrigen wird. In den 
ersten fünf Tafeln läuft n von 1 bis 60 Es sind auf- 
gezeichnet in TafelI: (1+i)r = ur (i= p/100), in 
Tafel II: derreziproke Wert davon 1/(1+i)r = vn, in 
Tafel III: a,)= (ur — 1):i, in TafelIV: ,];= (l— vr):i, 
in Tafel V: der reziproke Wert davon 1/a,| = Pz|-+ i, 
in Tafel VI: u112, il und 'ull2, in Tafel VII: 
j9)=12 (ur 2_1), jy—4 (ut -1) und jg=2(wrrR- 1), 
in Tafel VIIT: ö: j(,5), %: 3) und ö: j(), und endlich 
in Tafel IX: das Reziproke dieser Werte. Den ein- 
zelnen Tafeln sind zur Erläuterung ihrer Anwendung 
Beispiele beigegeben. Drei weitere zusammen- 
gesetztere Anwendungsbeispiele sind ausführlich am 
Schluß behandelt. Ausstattung und Druck dieser 
Tafeln ist ausgezeichnet. 


Dresden. 


Willers, 


Centre Belge de Recherches Mathematiques, Col- 
loque de G&omötrie algebrique. Tenn 
ä Lieges les 19, 20 et 21 Decembre 1949. 195 S. Liege 
et Paris 1950. George Thone, Masson & Cie, Preis 
200 fr. belg., 1400 fr. frangais. 


Bericht über ein internationales Colloguium vom 
19.—21. Dez. 1949 in Lüttich. Es hatte den Zweck 
einer Aussprache der Vertreter der verschiedenen Rich- 
tungen und Methoden in der algebraischen Geometrie. 
Im einleitenden Vortrage sprach F. Severi (Rom) 
über die italienischen Methoden und Probleme. Er 
untersucht dabei sehr tiefgründig die Legende von der 
fehlenden Strenge der italienischen algebraischen Geo- 
metrie; seine Gedanken gehen in ihrer Bedeutung 
weit übsr den Rahmen seines engeren Fachgebietes 
hinaus. Er unterscheidet zwischen substantieller und 
formeller Strenge und demonstriert, was er im Sinn 
hat, an zahlreichen Einzelproblemen. 


In den Vorträgen von Dubreil, v.d. Waer- 
den, Samuel kam dann die idealtheoretische 
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Richtung der algebraischen Geometrie zu Wort. 
Chätelet und Libois sprachen über die enge 
Verbindung und gegenseitige Befruchtung von Geo- 
metrie und Algebra im allgemeinen, B. Segre über 
die Beziehungen zu Zahlentheorie (diophantische Glei- 
chungen, Bureau über Anwendungen auf die 
Theorie der partiellen _Differentialgleichungen. 
Godeaux gab Anwendungen seiner Theorie der 
zyklischen Involutionen auf einer algebraischen 
Fläche. 

Das Collogquium und der vorliegende Bericht ist 


eine bedeutsame Tat der Verständigung und der Syn- 
these der verschiedenen algebro-geometrischen Schu- 


len. Dem fernerstehenden Mathematiker vermittelt 


er einen lehrreichen Einblick in Geist, Methoden und 
Fragen der gegenwärtigen algebraischen Geometrie. 


Dresden. Ott-Heinrich Keller. 


Dr.6.Bol,Elementeder Analytischen 
Geometrie. I. Gerade Linie und Ebene, 
Kreise und Kugeln, Kegelschnitte, alge- 
braische Kurven. 232 S. mit 103 Abb. Göttingen 
1948, II. Lineare, insbes. affine Transforma- 
tionen, Determinanten, Matrizen, Gruppen. 
156 S. mit 23 Abb. Göttingen 1949. Vandenhoeck u. 
Ruprecht. Preis kart. 14,— bzw. 8,80 DM, geb. 
16,50 bzw. 10,50 DM. 


Das Buch ist aus Freiburger Vorlesungen entstan- 
den. Der Stoff ist so ausgewählt und angeordnet, wie 
es den Bedürfnissen eines Studierenden der Mathe- 
matik oder Physik des ersten Semesters entspricht. 
Die Methoden sind bewußt elementar und ungekün- 
stelt gehalten. Verf. führt zunächst den Vektorbegrift 
ein und baut damit parallel die Geometrie der Ebene 
und des Raumes auf. Er geht überall vom Anschau- 
lichen aus und entwickelt allmählich die übergeord- 
neten Begriffe, wenn sich das Bedürfnis dazu ein- 
gestellt hat, und auch nur in der gerade nötigen All- 
gemeinheit und unter Verdeutlichung durch geschickt 
ausgewählte Beispiele. Durch Strenge und Gründlich- 
keit der Begriffsbildungen, durch Hinweise an allen 
passenden Stellen ist dabei der spätere Ausbau der al- 
gebraischen und axiomatischen Methoden vorbereitet. 
Jedem Paragraphen sind zahlreiche Übungsbeispiele 
beigegeben. Sehr sorgfältig und ausführlich sind die 
historischen Anmerkungen. Einem durchschnittlich 
begabten Anfänger wird das Studium des Buches keine 
besonderen Schwierigkeiten bereiten, und er wird nach 
gründlichem Durcharbeiten bemerken, welch tragfähige 
Grundlage er gelegt und welch sichere Beherrschung 
der Methoden er erworben hat. 


Ott-Heinrich Keller. 


Heyde, H. (Prof. a. d. Humbold-Universität Berlin): 
Mechanik für Ingenieure. BandI: Statik, 
Dynamik. (Teubners Fachbücher für Maschinenbau 
und Elektrotechnik.) X + 322 S. mit 317 Abb. Leipzig 
1950. B. G. Teubner Verlagsgesellschaft. Preis geb. 
13,40 DM. 


Die starke Nachfrage nach guter Literatur über 
technische Mechanik, die sich nach Kriegsende be- 
merkbar machte, hat inzwischen durch verschiedene 
hier besprochene Neuauflagen und Neuerscheinungen 
wesentlich nachgelassen. Das vorliegende Buch gehört 
nicht zur wissenschaftlich-formalistischen, sondern zur 
rein ingenieurmäßigen Mechanik-Literatur. Es kam 
dem Verf. offenbar in erster Linie auf die Darstellung 
der einfachen Anwendungsformeln nach Art eines 
Taschenbuches an, wobei man den Verzicht auf ge- 
wisse theoretische Grundgedanken nachteilig empfin- 
det. Während die Darstellung des ersten Teiles (Sta- 
tik) noch Einfachheit und Klarheit aufweist, macht 
sich im zweiten Teil (Dynamik) das Fehlen analyti- 


scher Bestandteile der Mechanik störend bemerkbar, 
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So werden die Schwingungsprobleme ohne Differential- 
gleichung als Beispiel der rein kinematischen Bezie- 
hungen des Bewegungsvorganges behandelt. Dadurch 
verliert der Leser das Gefühl für die Kausalität und 
die Reihenfolge der beim deduktiven Denken zur An- 
wendung kommenden Gesetze. Die angenommene 
Wegzeitfunktion und ihre Ableitung werden gewisser- 
maßena priorials ,, Bewegungsgesetz‘‘ und ‚‚Geschwin- 
digkeitsgesetz‘‘ deklariert, während die dynamische 
Grundgleichung ganz in den Hintergrund tritt (das 


D’Alembertsche Prinzip bleibt ignoriert). Durch eine 
solche Darstellungsart wird vielen Lesern der Weg zur 
höheren Mechanik von vornherein verschüttet; denn 
sie werden nicht an die in der Mechanik erforderliche 
besondere Art des deduktiven Denkens gewöhnt; Die 
Erziehung dazu muß aber als eine der wesentlichsten 
Aufgaben aller Mechanik-Lehrbücher angesehen wer- 
den. 


Dresden. H. Neuber. 


EINGEGANGENE BÜCHER 


Bei der Schriftleitung sind folgende Bücher eingegangen (ausführliche Besprechung bleibt vorbehalten). 


Dr. phil. A. Duschek (o. Prof. a.d. techn. Hochschule 
TH. Wien, Vorlesungen über höhere 
Mathematik, 2. Band, VI-+ 386 S. mit 125 Abb. 
Wien 1950. Springer-Verlag. Geh. 26,— DM., geb. 
29,— DM. 


Dr.-Ing. R. €. Oldenbourg und Dr.-Ing. H. Sartorius, 
Dynamik selbsttätiger Regelungen. 
Band I. Allgem. und mathem. Grundlagen. Stetige 
und unstetige Regelungen. Nichtlinearitäten. 2. Aufl. 
2588. mit 112Abb. und 1 Tafel. München 1951, 
Verlag R. Oldenbourg. Preis geb. 26,— DM. 


H. W.v. Gusrard (Direktor des Statist. Amtes der 
Stadt Düsseldorf), Über Zahlen und’Dia- 
gramme zur Messung des Verkehrs 
(Wissensch. Veröffentlichungen und Tagungsberichte 
des Minist. für Wirtschaft und Verkehr des Landes 
Nordrhein-Westfalen, Heft 12.) 378. mit 4Abb. 
Düsseldorf 1950. Droste-Verlag. Preis 2,50 DM. 


A. Jante (o. Prof. a. d. Techn. Hochschule Dresden), 
Leitfaden der technischen Thermodynamik. 
(Teubners Fachbücher für Maschinenbau und Elektro- 
technik.) VI+146 S. m. 153 Abb. Leipzig 1950. 
Teubners Verlagsgesellschaft. Preis geb. 9,60 DM. 


P. Crantz, Sphärische Trigonometrie. Bearb. 
von Dr. M. Hauptmann. (Mathematisch-physika- 
lische Bibliothek, Reihe II, 7.) 4. Aufl., 112 S. m. 76 Fig. 
und 1 Nomogramm. Leipzig 1950. B. G. Teubners 
Verlagsgesellschaft. Preis kart. 4,20 DM. 


Dr. phil. habil. W. Gebelein (Landshut) und 
Dr. med. H.-J. Heite (Münster), Statistische 
Urteilsbildung. Erläutert an Beispielen aus 
Medizin und Biologie. Mit einem Geleitwort von 
Prof. Dr. med. C. Moncorps (Münster). 192 S. 
mit 50 Abb. und 20 Beispielen. Berlin- Göttingen- 
Heidelberg 1951. Springer-Verlag. Preis: brosch. 
13,60 DM. 


Die besprochenen und angezeigten Bücher sind durch den Buchhandel zu beziehen, 


iR NACHRICHTEN 


Gießen: Profesor Rolf Nevanlinna, Hel- 
singfors und Zürich, wurde das Ehrendoktorat der 
Naturwissenschaften von der Naturwiss. Fakultät der 
Justus-Liebig-Hochschule Gießen verliehen. 

Ende Februar 1951 fand im Mathematischen Insti- 
tut der Justus-Liebig-Hochschule Gießen eine funk- 
tionentheoretische Arbeitstagung statt, bei der vor 
allem jüngeren Funktionentheoretikern aus Karls- 


ruhe, Marburg, Münster, Zürich und Gießen Gelegen- 
heit zu gründlichem Kennenlernen neuester Arbeiten 
und zum Bericht über eigene Forschungen gegeben 
war. 


Hannover: Am 11. Febr. 1951 verstarb der 
a. o. Professor an der T, H. Hannover, Dr. Conrad 
Ludwig. 


ZUSCHRIFIEN AN DEN HERAUSGEBER 


Zu W. Schmidt, Untersuchungen des für 
den schiefen Stoß elastischer Kugeln gültigen 
Reflektionsaxioms und einige Folgerungen 
daraus. Z. angew. Math. Mech. 30 (1950), S. 182 
bis 184. 

Herr W. Schmidt hat in den kleinen Mitteilungen 
dieser Zeitschrift eine Note mit dem Titel veröffent- 
lieht: Untersuchungen des für den schiefen Stoß 
elastischer Kugeln gültigen Reflexionsaxioms und 
einige Folgerungen daraus. Er behauptet mit Recht, 
daß die oft angenommene Gleichheit von Einfalls- und 
Ausfallswinkel im allgemeinen falsch ist. Das steht 
auch schon in der Lösung der Aufgabe, die man in 
meiner ‚‚Elementaren Mechanik“ unter Nr. 298 in 
852, S.448 findet und zwar wesentlich allgemeiner, 
nämlich unter Einbeziehung der Gleitreibung, die 


Herr Schmidt ausschließen will. 


Doch unterscheiden sich unsere Betrachtungen 
weiterhin in einem wesentlichen Punkt. Herr Schmidt 


hat noch eine zweite Hypothese für die tangentiale 
Komponente des Stoßes und somit eine zweite Kon- 
stantegmit 0<g?<l. Das würde eine Gleichung zu 
viel geben, wenn nicht Herr Schmidt von der Bedin- 
gung absähe, daß die tangentiale Geschwindigkeit am 
Ende null sein muß, wie es der ausdrücklich von ihm 
formulierten Voraussetzungen entspricht, daß nur 
Haftreibung, keine Gleitreibung vorhanden sein soll. 

Vergleicht man unsere Resultate, so stimmen sie 
überein, wenn man q= 0 setzt. 

Aber meine Theorie ist deshalb nicht ein Sonderfall 
von der des Herrn Schmidt. Ich mußte diese zu- 
nächst für falsch halten, da sie etwas zu wenig und 
etwas zu viel enthält. Aber ich finde jetzt, daß sie 
doch ihren Sinn hat, sie behandelt nur ein anderes 
Problem. 

Wenn man den Stoß einer Kugel gegen eine feste 


Wand behandelt, betrachtet man diese in der Regel, 


so wie ich, als makroskopisch eben. Dann hat man 
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bei Annahme von Haftreibung die Bedingung, daß 
die tangentiale Geschwindigkeit zum Schluß null'sei, 
d.h., daß in der Bezeichnung des Herrn Schmidt 


sein muß. Dann aber darf nur eine Stoßhypothese 
gemacht werden, etwa ebenso wie bei Herrn Schmidt 


BE (8), 


r w 
vosinB—— —- =0(.. 
C 


vcos B=kucos A 
Darin und in der Impulsgleichung 
vısinBtwe=usnA ...:. 


stimmen wir überein, 
(2), (8) und (a) genügen, um »v, B und w zu be- 
rechnen. Es ergeben sich 


V—IU Vi: cos? A mi sin? A, 
1 
tg B Tanne 
Se T 
w=usin A Iire' 


Wegen der Bedingung der Haftreibung kommt noch 
eine Ungleichheit hinzu, die wir nicht hinschreiben 
wollen. 

Statt (a) hat Herr Schmidt die zweite Stoßhypothese 


. (b), 


ohne sie explizit hinzuschreiben. (a) ist zwar mathe- 
matisch ein Sonderfall von (b), hat aber einen ganz 
anderen Charakter, Sie ist kinematischer Natur, 
lediglich der Voraussetzung des Haftens entspringend, 
während (b) eine Stoßhypothese ist. 

Das braucht nun doch nicht falsch zu sein, wenn 
man sich die Wand und vielleicht auch die Kugel als 
so rauh vorstellt, daß sie sich beim Stoß ineinander 
verhaken, man also mit der Annahme einer einheit- 
lichen wohlbestimmten Berührungsfläche nicht mehr 
auskommt und so veranlaßt sein kann, der Un- 
genauigkeit des Problems entsprechend eine zweite 
Stoßhypothese einzuführen, den tangentialen Stoß 
also gesondert zu betrachten. Man hat gewissermaßen 
zwei Wände. Setzt man nun g nicht gleich null, so 
muß die Annahme der Haftreibung fallen, denn die 
Kugel verläßt die Wand mit einer von null verschie- 
denen Geschwindigkeit, welche der rechten Seite 
von (b) entspricht. Wie sie das freilich bei der Ver- 
zahnung fertig bekommt, ist noch eine Frage, Aber 
immerhin mag für dieses undurchsichtige Problem die 


} w ; 
vesinB—— =—qusinA.. 
6 


Betrachtung des Herrn Schmidt Wert besitzen: zur‘ 


Anpassung an die Wirklichkeit steht ihm ja noch 
das q frei. 

Im Grunde behandelt also Herr Schmidt ein ganz 
anderes Problem als das, was man sich unter dem 
Stoß gegen eine feste Wand vorzustellen pflegt. 
Offenbar versteht auch Herr Schmidt unter Haft- 
reibung etwas anderes als man es in der Regel tut. 
Er meint wohl, daß die tangentiale Geschwindigkeit 
das Zeichen wechselt und also auch einen Augenblick 
null ist, 


Landshut (Bayern). G. Hamel. 


Erwiderung 

Die freundliche Zuschrift von Herrn Prof. Dr, 
G. Ha mel gibt mir Gelegenheit zu notwendigen Er- 
gänzungen meines Aufsatzes. Erst nach dessen Ver- 
öffentlichung wurden mir Schriftquellen zugänglich, 
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 sammengesetzt und miteinander vermischt, daß sie 


„verständlich sind ebenso wie g = — 1 (reibungsfreies, 


die im Gegensatz zu den apodiktischen, aber falschen‘ 
Aussagen über das Reflexionsaxiom beim schiefen 
Stoß stehen, bei denen aber trotzdem für den Eh 
meinen Fall Ungleichungen nicht entbehrt werden 
können und die zudem durchweg die Realität dr 
tangential wirkenden Federkraft außer Ansatz lassen. 

Daß ich entgegen der Ansicht von Herrn Prof. Dr. ) 
G.Hamelin meinen Betrachtungen die Gleitreibung 
nicht ausschließe, folgt zwingend aus dem Schlußsatz 
meiner Abhandlung, in dem ich mich der Fiktion 
eines rein körperlichen Stoßvorganges bei der Reflexion 
des Lichtes bediene. Bei dieser Fiktion wird also 
qg=— 1. Es handelt sich hierbei, wie betont, um 
den Grenzfall völlig reibungsfreien Gleitens — selbst- 
verständlich also auch völlig haftfreien Gleitens —, 
wie er vielfach fälschlich noch heute im Fachschrift- 
tum zur Erklärung des Reflexionsaxioms beim schiefen 
Stoß einer Kugel gegen eine Wand benutzt wird.‘ 
Unter Haftreibung verstehe ich mit Herrn Prof. Dr. 
G.H am elgenau das gleiche, wie z. B.in Sommerfeld 
„Mechanik“ 1944, S. 76ff. erläutert. Gerade von hier 
aus wird aber der Unterschied der beiderseitigen Be- 
trachtungen erkennbar. Wie in meinem Aufsatz, 
wenn auch kurz und unvollkommen, erläutert, super- 
poniereich die Kraft der Haftreibung sowohl wie 
auch die der Gleitreibung der tangential wirkenden 
Federkraft, die sich entsprechend den zu Beginn 
meines Aufsatzes erläuterten Versuchsergebnissen b- 
merkbar macht. Zusammenfassend sind die beim Ee 
schiefen Stoß in ihrer Gesamtheit wirkenden Kräfte, 
insbesondere also Haftreibung, gleitende Reibung und 
Federkraft, ihrem Ursprung nach so verwickelt zu- 


noch weniger exakter Analyse zugänglich sind wie 
diejenigen beim Zentralstoß. Um so empfehlens- 
werter ist in Analogie zur Newtonschen Stoßzahl .k, 
die sich ausschließlich auf die Normalkomponente des 
Stoßes bezieht, die Einführung der in meinem Aufsatz 
vorgeschlagenen Stoßzahl q für die Tangential- 
komponente. Übrigens schließt k offensichtlich auch 
Reibungsvorgänge in sich ein; diese Reibungskräfte _ 
entstehen beim Äneinandergleiten der sich gegenseitig 
verformenden Stoßflächen, heben sich aber in diesem a 
Fall in ihrer richtungändernden dynamischen 2 
Wirkung sicher weitgehend gegeneinander auf. Selbst- 


daher auch federkraftfreies Aneinandergleiten) auch 
die Werte g = + 1 (verlustfreies, gleitfreies Haften in 
Verbindung mit verlustfreier tangentialer Federkraft) 
und gq = 0 (verlustfreies, gleitfreies Haften ohne Wirk- 
samwerden der tangentialen Federkraft) nur Ideal- 
fälle. In Wirklichkeitistt +1>g>— Inurempirisch 
zu ermitteln, genau wie k; sicherlich ist g auch ab- 
hängig von d, uw und A, genau wie erwiesenermaßen 
auch k von u abhängigist. qläßt sich genau wie % fall- 
weise experimentell einwandfrei ermitteln. Für die 
mathematische Behandlung der Vorgänge sehe ich 
nach wie vor in der eleganten Poinsot-Methode (Ein- 
führung der reduzierten Massen) das einfachste, an- 
schaulichste und aufschlußreichste Verfahren. Wenn 
man die Vorgänge beim schiefen Stoß genauer unter- ” 
suchen will, bietet die vorgeschlagene Tangentialstoß- 
zahl q m. E. eine ebenso wertvolle Stütze wie k, klärt 
die Zusammenhänge und ist eine notwendige Ergän- 
zung der analytischen Hilfsmittel für die heute im: - 
noch recht problematischen Vorgänge beim sch 
Stoß, allerdings keine hinreichende, sofern, wie 
Aufsatz angedeutet, die Schwerkraft und die Stoß- 
wirkung der Gezeiten in die Betracht mit ein- 
bezogen wird, 28 

Rheydt. 


